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Abstract. Colored tensor models have recently burst onto the scene as a promising con- 
ceptual and computational tool in the investigation of problems of random geometry in 
dimension three and higher. We present a snapshot of the cutting edge in this rapidly 
expanding research field. Colored tensor models have been shown to share many of the 
properties of their direct ancestor, matrix models, which encode a theory of fluctuating two- 
dimensional surfaces. These features include the possession of Feynman graphs encoding 
topological spaces, a 1/N expansion of graph amplitudes, embedded matrix models inside 
the tensor structure, a resumable leading order with critical behavior and a continuum large 
volume limit, Schwinger-Dyson equations satisfying a Lie algebra (akin to the Virasoro al- 
gebra in two dimensions), non-trivial classical solutions and so on. In this review, we give 
a detailed introduction of colored tensor models and pointers to current and future research 
directions. 
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1 Introduction 

Many problems one comes across in mathematics and theoretical physics, after one strips away 
the interpretational dressing, reveal themselves to be enumerative in nature. Simply speaking, 
one counts. Meanwhile, the advent of quantum theory revealed the evermore growing importance 
of probabilistic processes in physics. Much ink has been spilled in the mathematical community 
pursuing the development of probability theory in terms of measure-theoretic concepts [99]. 
Far from being divorced, these combinatorial and probabilistic concepts appear to facilitate 
a harmonious co-existence. In fact, it is rather the case that they are becoming increasingly 
intertwined. Historically, the model processes of probability theory often involved an explicit 
combinatorial mechanism, for instance, counting the ways of placing balls of varied colors into 
urns [27]. These seemingly antiquarian methods had direct physical implications relating to 
contemporary voting procedures and models of gases. One could perhaps view the theory of 
particles undergoing Brownian motion as the succinct culmination of these efforts. It also serves 
to illustrate these probabilistic approaches in a more modern setting: summing the (random) 
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walks of the particles with a probability measure captures the diffusive nature of this physical 
process. 

This marriage provides powerful tools to investigate a multitude of theories, nowhere more 
prevalently than in quantum field theory. In this context, one has a priori three approaches 
within which to pose a given problem: Schrodinger's analytic approach, Heisenberg's algebraic 
approach or Feynman's combinatorial approach. While certain conceptual problems may favor 
solution by analytic or algebraic methods, it is undoubtedly true that the combinatorial for- 
mulation has proved to be the preeminent instrument for extracting data. One of its strengths 
as a computational tool stems from the fact that many applications just involve the ability 
to perform (generalized) Gaussian integrals. Furthermore, one need look no further than lat- 
tice QCD to find a scenario very well adapted to numerical methods in the non-perturbative 
regime. 

Meanwhile, the combinatorial approach also has many attractive features from a conceptual 
viewpoint: it implements the symmetries of the theory directly; it allows one to incorporate 
constraints in a simple fashion and thus isolates relevant dynamical variables. In fact, many 
thermodynamical quantities (free energies, critical exponents) and observables (Wilson loop 
observables) have their most natural definition in this language. 

In this review, we shall discuss a class of quantum field theories known as colored tensor 
models. Since we develop the theory in quite a linear fashion, we feel it is appropriate to 
motivate each section individually. 

Gaussian measures and Feynman graphs - with what kind of theory are we dealing? The 
combinatorial formulation of quantum theory has exercised such an appealing allure that here, as 
in many other texts, we shall define a quantum field theory as a probability measure for random 
functions. In the spirit of the physics literature, we shall call the random functions 'fields'. We 
start from rather general principles and introduce the theory of probability measures adapted to 
our quantum field theoretic setting, that is Gaussian measures and their perturbation. Having 
said that, we do not stray into the more arid regions of measure theory involving fi-fields, 
Borel sets and the like (for that, we invite the reader to consult any thorough textbook on 
probability theory, for example [76]). Throughout this section, we supply ample simplified 
examples, relevant for out subsequent study of colored tensor models. 

The physical observables are the correlations of the probability measure. Whenever the 
probability measure is a perturbed Gaussian measure, the observables of a quantum field theory 
are computed via an expansion in Feynman graphs. The perturbation is encoded in an action, 
a polynomial functional of the random fields. The graphs, which are combinatorial objects, 
come equipped with combinatorial weights and amplitudes fixed by the details of the action: 
the measure and the vertex kernels of the perturbation. The graphs comprise of vertices and 
edges, hence topologically, they are 1-complexes. 

To conclude this section, we introduce the generic form of {D + l)-colored quantum field 
theory models. They are probability measures for a collection of random fields indexed by 
colors. In turn, this leads to graphs in which the edges possess a color index {{D + l)-colored 
graphs) . 

A brief review of (colored) matrix models - where can we look for hints? In this section, we 
present a quick review of matrix models, emphasizing the aspects that we shall later generalize 
to colored tensor models. In the case of matrix models, the coloring does not play an important 
role and most of the notions we present for the colored models are in fact relevant to non-colored 
ones. However, in order to prepare the reader for the higher-dimensional case, in which the colors 
play a crucial role, we highlight the colored case also for matrices. 

Topology of colored graphs - what are the properties of our Feynman graphs? We examine 
in detail the characteristic properties of the (D + l)-colored graphs. The nested structure of 
subgraphs indexed by colors endows the graphs with a Z'-dimensional cellular complex structure. 
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Most importantly, each cell complex associated to a graph is a simplicial pseudo-manifold, 
consequently the colored models are statistical theories of random i^-dimensional topological 
spaces. We present the colored cellular homology and homotopy for our graphs and clarify the 
structure of their boundaries, relevant for the analysis of observables of the colored models. 

Subsequently, we proceed to the core of this section. From the lessons one learns in the study 
of matrix models, it is of utmost importance that one identifies combinatorial and topological 
quantities associated to graphs. For matrix models, this is the genus of the (ribbon) Feynman 
graph and of its dual Riemann surface. Something very similar happens in higher dimensions for 
the tensor models we study. For any colored graph there is a precisely defined set of Riemann 
surfaces embedded in the graph, known as jackets. Each such jacket corresponds to a ribbon 
graph embedded in the {D + l)-colored graph. Each jacket has then a genus, and one defines 
the degree of the (D + l)-colored graph as the sum of all the genera of its jackets. 

The picture, however, in higher dimensions is not as rosy as it might seem so far. The degree 
is not a topological invariant and the study of the full perturbative expansion of tensor models 
requires more care. It transpires that there are a set of graph manipulations which preserve 
certain combinatorial properties: these are the A:-dipole creation/annihilation moves which were 
already known in the theory of manifold crystallization [64, 104]. They allow one to define 
the notion of combinatorial equivalence and classify graphs into equivalence classes. Moreover, 
a subset of fc-dipole moves encodes homeomorphisms and thus preserve the topology of the cell 
complex associated to the graph. This subset leads to a classification of graphs taking into 
account not only the combinatorics but also the topology. 

Our interest in the fine structure of jackets does not end there. It appears that one can 
partition the set of jackets into subsets which individually capture the degree of the graph. This 
partition reappears when dealing with embedded matrix model regimes and classical solutions. 

Tensor models. To a certain extent, the preceding three sections were mathematical prelimi- 
naries. In particular, no reference has so far been made to any specific higher-dimensional tensor 
model. We shall now deal with the first such specific example. A statistical theory of random 
topological spaces requires no other ingredient than the colors. On the contrary the weights 
(amplitudes) of the graphs are strongly dependent on the details of the model: the number of 
arguments of the fields, the connectivity of the arguments at the vertices, the Gaussian measure 
with respect to which we perturb, an so on. We will concentrate our study on the simplest 
model one can consider, the independent identically distributed (i.i.d.) colored tensor model. 

Before we go any further, let us spare a moment to place colored tensor models within the 
larger scheme of things. Often the random fields in quantum field theory are defined on a group 
and one can formulate the same quantum field theory in terms of its Fourier modes. The random 
field thus becomes a random tensor. As M'^ is a group, most of the quantum field theories one 
encounters in flat space time (quantum electrodynamics, non-Abelian gauge theories, etc.) can 
be considered, in the broadest sense, tensor models. In opposition to this, we will call tensor 
models in this review those quantum field theories whose graphs encode a topological space of 
dimension at least two. The simplest tensor models are thus the random matrix models. We 
will use synonymously the name group field theories for these models^. 

Much of the analysis we conduct in this paper is the generalization in higher dimensions of the 
one undertaken in matrix models (which are probability measures of random rank two tensors) . 
Random matrices first appeared in physics as an attempt to understand the statistical behavior 
of slow neutron resonances. Later, they were used in efforts to characterize chaotic systems, 
describe elastodynamic properties of structural materials, typify the conductivity attributes 
of disordered metals, study the theory of strong interactions [128], detail aspects of various 
putative physical systems such as two-dimensional quantum gravity (both pure or coupled to 



'^Less inclusive authors prefer to use the name group field theory when referring to specific subclasses of such 
models. 
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matter) [39, 41, 46, 53, 93, 94, 98], conformal filed theory [47, 57, 60], string theory [43, 59, 77], 
illustrate the distribution of the values of the Riemann zeta function on the critical line, count 
certain knots and links and the list keeps growing. In all these applications, the power of the 
approach resides in the control one has on the statistical properties of large matrices, of size N. 
This level of control stems from the so called expansion', which is a (much better behaved) 
alternative to the usual small coupling expansion in quantum field theory dominated by graphs 
corresponding to the simplest (spherical) topology [42, 46, 93]. 

Like matrix models, colored tensor models come equipped with a large parameter N, the 
size of the tensors. It is of utmost importance that we identify in this case also an expansion 
in and understand it as thoroughly as possible order by order. It is only then that we have 
any chance of comprehending the statistical behavior of random topological spaces in higher 
dimensions. 

The material discussed in this rest of this section is dedicated to presenting the 1/A^ expansion 
of colored tensor models. In fact in higher dimensions we must distinguish between two distinct 
expansions, one taking into account only the amplitude of the graphs and one including also 
their topology. To this end we show that the power counting in of the amplitudes of the 
i.i.d. model is given by their degree. We then make use of the combinatorial and topological 
moves to define the notion of core graph. Roughly speaking, these are the 'simplest' graphs 
(according to specific criteria) of a certain degree in a combinatorial/topological equivalence 
class. The power of these core graphs is that one can use them to order the 1 /N expansion. We 
give a rundown of several redundancies in these expansions and catalog the lowest order examples 
of core graphs. Critically, the leading order graphs correspond to the simplest, spherical topology, 
in any dimension. However, not all graphs of spherical topology, but only a very specific subclass, 
contribute at leading order. 

Embedded matrix models. By the time we get to this point of the review, we shall have heard 
much about the similarities and differences between matrix and tensor models. There are certain 
algebraic and analytic tools we would like very much to generalize from the matrix to the tensor 
scenario vis a vis their exact solution [53] (at least in a certain regime). We take a first step in 
this direction here. We utilize the concept of PJ-factorization to partition the degrees of freedom 
of the tensor model. Rather than viewing the tensor model as a simply generating weighted D- 
dimensional cellular complexes, we can identify already at the level of the action, the generators 
of the embedded jackets. These generators correspond to matrix models embedded inside the 
tensor model. This provides a launch pad to applying matrix model techniques directly to the 
certain tensor probability measures. 

Critical behavior. In this section we investigate in detail the dominant contribution in the 
large N limit of colored tensor models. We provide a purely combinatorial characterization of the 
'melon graphs' contributing to the leading order generalizing the planar graphs [42] to arbitrary 
dimensions. By constructing an explicit map between graphs and colored rooted (D + l)-ary 
trees, we present several ways to resum the series. We show that this resumed leading or- 
der exhibits a critical behavior and undergoes an analytically controlled phase transition from 
a discrete to a continuum theory. We subsequently present an interpretation of the i.i.d. col- 
ored tensor model in terms of dynamical triangulations. In this more geometric perspective 
the leading order continuum phase shares the critical behavior of the branched polymer phase 
of dynamical triangulations. However, the average spectral and Hausdorff dimensions of the 
melonic graphs ensemble have not been computed and one can not yet conclude on the precise 
relationship between melonic graphs and branched polymers. 

Bubble equations. A key objective of any quantum field theory is to identify the quantum 
equations of motion, that is, the equations satisfied by the correlation functions. These are the 
Schwinger-Dyson equations and contain all the information pertaining to the quantum dynamics. 
Moreover, they can be re-interpreted as operators on the space of observables. These operators 
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form an algebra, just one representation of which is given by the quantum field theory correlation 
functions. This opens the door to finding other faithful representations of the algebra, which 
will describe essentially the same theory, although their aesthetic presentation might emphasize 
different aspects. It is in this manner that one maps between matrix models and Liouville 
gravity in two dimensions gravity since the Schwinger-Dyson (also called 'loop') equations form 
a representation of a sub-algebra of the Virasoro algebra [13, 71, 108]. 

Classical solutions. Any analysis of a field theory would be incomplete if one did not attempt 
to uncover pertinent information about its classical regime. For some field theories, such as 
(quantum) electrodynamics, the classical action describes directly rich physical phenomena. For 
others, such as quantum chromodynamics, the classical action holds the key to unlocking non- 
perturbative information. For example, instantonic solutions furnish genuinely non-perturbative 
gauge field configurations which display a host of geometrical, topological and quantum effects 
with fundamental impact for the ground state and spectrum of non-Abelian gauge theories. 
While we do not undertake an exhaustive study of the classical theory of colored tensor models, 
we shall analyze a class of ansatze for solutions of the equations of motion. 

Subsequently, we shall perturb the theory in various ways about this non-trivial solution and 
investigate the resulting quantum behavior. 

Extended discussion and conclusion. We have not attempted here to present a complete 
literature review. In particular, we have chosen not to include many interesting developments 
that have taken place in specific tensor models, but which did not fit into the flow of the main 
review. For that reason, we have expanded the discussion section so that we could at least 
briefly coordinate these results with what we have presented here. 

Much remains to be studied about colored tensor models, its relation with higher-dimensional 
conformal field theory, statistical physics and quantum gravity, analysis of sub-dominant con- 
tributions, tighter control over the topological features of the graphs, just to name the first few 
that spring to mind. We tender a forward-looking review to conclude. 



2 Gaussian measures and Feynman graphs 

Ultimately, we shall be interested in analyzing a class of quantum field theories known as colored 
tensor models. To this end, we first define what we mean by a quantum field theory and present 
the tools required to investigate its properties. 

In the first part of this section, we review the relevant aspects of perturbed Gaussian prob- 
ability measures, along with the Feynman graphs occurring in the perturbative expansion of 
the free energy. In the second part, we introduce the colored measures on which colored tensor 
models rely. Most of the material we present is standard. 

Let X be a nonempty set and denote its elements hy n £ X. We consider (/> : X — )• M (or C, 
or G), a real, complex or Grassmann- valued random function defined on X.^ H is a complex 
(Grassmann) function, we denote (j) its complex conjugate (its involution). By ^iX), we mean 
the set of all functions on X and by 0„ = (pin) the value of (p at the point n. 

Definition 2.1. A quantum field theory is a probability measure — dv{(f)) for a random function 
(field) (j), where Z is the normalization factor. 

In the sequel, we shall identify a probability measure through its partition function Z and 
correlations: 

^ = (l) = y" du{ci)), {(j)n^ ...(f)^^) = j dv{4>) (j)ni--- 4>np- 

^We consider only Grassmann algebras endowed with an anti-involution <p ^ (f> with 0x 



= -x0, <t> = -4>- 
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The logarithm of the partition function is called the free energy, and is denoted F = \nZ. 
Importantly, generic correlations may be expressed in terms of a subset known as connected 
correlations (and signified by {■)c)- 



where P denotes the partitions of the set {1, . . . ,p} into non-empty subsets S. 

Example 2.1. The set X has a unique element X = {0}. The random function (po is then 
a random variable. A skewed coin with (real-numbered) outcomes a and b, for instance, is the 
quantum field theory of the pure point probability measure: 



Example 2.2. The set X has a finite number of elements X = {1, . . . ,N}. A quantum field 
theory is encoded in a probability density P: 



where dcpn is the usual Lebesgue measure on M. 

Of course when the cardinality of X becomes infinite, the definition of a probability measure 
is a thorny issue. However, normalized Gaussian measures, can be defined simply also in this 
case. Suppose that X is a measurable space, and label the Lebesgue integral over X by J2n- 
A linear operator C : ^{X) — >■ ^{X) is identified by its kernel C„^„/: 



For complex fields, we will sometimes denote their variables by n. Both n and n belong to X, 
that is, for the elements, the bar does not denote complex conjugation; it is just a bookkeeping 
device used to track the indices belonging to a complex-conjugated field. 

Example 2.3. Consider X the set: 

X = {n\n = (m, . . .yUo), nfc = 1, . . .,N}. 

The Lebesgue measure on X is a pure point measure. A random function on X is a random 
tensor with D indices ^fj. The kernel of an operator C is a A^^ X A^^ matrix Cnn'- 

Definition 2.2. A normalized Gaussian probability measure of covariance C is a measure whose 
only non-zero correlations are: 

• for a real field ^ : X ^ M, 



p seP kes 



diy{4'o) = d(po{pSia- (po) + qS{b- (po)), Z ={p + q), 



{^^)=pa^ + qb\ 






where V denotes all the pairings (that is, distinct partitions of the set {!,..., 2p} into 
2-element subsets), and {i,j) denotes a pair. 
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• for a complex field (f)^(j) : X ^ 

\fc=l / TT k 

where tt is a permutation of p elements. 

• for Grassmann fields <j), (j) : X ^ G, 

V 

k=l vr k 

where e(7r) is the signature of the permutation vr. 

The covariance C needs neither be an invertible nor a hermitian operator. 

Example 2.4. \i X = {!,..., A^}, then the covariance is a matrix. Suppose the covariance 
is invertible and denote its inverse by C~} . The normalized Gaussian probability measure of 
covariance C can be written: 

• for a real field : X — t- M, 

N 
n.Ji' — 1 



n=l ^ 



for a complex field (p : X ^ C 

N 

dM'A>) = det(C)n^^e 

n= 

for a Grassmann field 6 : X 



N 



n=l 



n=l 



n,n— 1 



where, in this final case, the ordering matters^. 
A polynomially perturbed Gaussian complex or Grassmann measure is: 

oo 
s,s=0 

where the repeated indices n and n are summed (note that s and s are independent). 
Example 2.5. Let X = {0} and (pQ a complex random variable. The D-ary tree measure is: 

duJ^s(^(t>oAo) = dfii{^o,^o)e-^^'^+9<l>oi^. 



^The integral over the Grassmann algebra is defined as f d(j) = d^, and the derivative is defined 'to the left' 
d^iM) = A, iff A 7^ (PB. 
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Example 2.6. The complex {or Grassmann) measure is: 



= dnc{(l>,4>)e~ 



The partition function and correlations of perturbed Gaussian measures are evaluated through 
Feynman graphs. These graphs are obtained by Taylor expanding with respect to the perturba- 
tion parameters (the coupling constants) and computing the Gaussian integrals. The correlations 
of t fields 6 and i fields d> can be written as: 



(0, 



ni Ynt 



Pa, 3=0 



n 

s,s 



Ps 



/=1 



/Ps,s 



(0 -(0 -(0 



where again all the repeated n, n indices are summed. The Gaussian integrals are zero unless 
t + Ylips s = i + Ylips ^- Denoting generically the n indices by n^, the n indices by n^, 
and TT a permutation of /c = t + __^q s elements: 



(0ni • • • 



nt Yni 



EE n 

Pa, 3=0 TT s,s 

Ps,s¥=0 



(i) (0 _(0 -(I) 
P= f! -^-^ n\'...ny;n\'...ny 



,7r(nfe)- 



Each term in the above sum can be represented as a graph. The kernels K are represented as 
(s + s)-valent vertices, having s half-lines (p and s half-lines (p. The insertions (pn and (pn are 
represented as external points. A permutation tt corresponds to a contraction scheme in which 
one connects all indices n with indices n. Every connection is represented as a line. The sum 
may be rewritten as a sum over Feynman graphs with ps^s vertices which are {s + s)-valent: 



>=E n 



-A 



ni...ne;ni...nj. 



A 



ni...ne;ni...ng 



(5)= n 

s,s 



Ps,s 



where A(Q) is the amplitude of the graph Q (again recall that repeated indices are summed). 
The class of graphs Q one sums over is fixed by the terms in the action. The free energy is the 
sum over connected graphs. Unless otherwise specified we always deal with connected graphs. 

Example 2.7. Consider the L>-ary tree measure. The partition function is: 

Z = [ d,xi{<PA)e''^+^^~^'' = V -^JV / di,i{4>,4>^^^'^''' 



p,g=0 
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and its 1-point correlation function is: 

p,q=0 q=0 V y / 

The graphs contributing to the connected 1-point function {4>)^ ~ 'E {^) connected graphs 
built as follows. They have a root vertex corresponding to (j). The root either connects to 
a 1-valent J- vertex (in which case, the graph consist in exactly one line) or it connects to the 
(j) field in a (D + l)-valent g-veitex. All the other fields of the (/-vertex are <j). In turn, each of 
them either connects to a J- vertex or to another (7- vertex. The graphs can not possess any loop 
line, since at every step in this iterative process, all as yet uncontracted fields are (f). Thus, the 
graphs contributing to the connected 1-point function of the D-ary tree measure are exactly the 
rooted D-ary trees. For further reference we note that: 

see [33]. 

Although one obtains a summable series in the case of the D-ary tree measure, in general the 
expansion in Feynman graphs is divergent. First of all, for numerous models, the amplitudes of 
individual graphs diverge. These are the well known ultra-violet divergences. In order to render 
the amplitudes finite, one introduces cutoffs and a renormalization procedure. A renormaliza- 
tion procedure relies on the splitting of normalized Gaussian measures. We shall not enter here 
in the details of this procedure (see [121] for an elementary introduction to the topic of renor- 
malization). What we stress is that the effective theory obtained following a renormalization 
procedure is unavoidably an effective theory for the IR degrees of freedom corresponding to 
large eigenvalues of the covariance. 

After the UV divergences are dealt with, one needs to address the issue of summability of 
the series. In generic models, one finds loop lines in the Feynman graphs. Due to the presence 
of these loop lines, the number of graphs grows faster than factorially (super exponentially) and 
the series has zero radius of convergence, being at most just Borel summable. 

The Feynman graphs of the simplest quantum field theories are made of vertices (0-cells) and 
lines (1-cells); hence, they are 1-dimensional cellular complexes. Tensor models are quantum 
field theories such that their graphs have a richer (combinatorial) topological structure, namely 
they also possess faces (2-cells). 

2.1 The {D + l)-colored models 

We come now to the definition of a colored model. As we will see in a later section the graphs 
of that colored model with D colors automatically have a D-dimensional cellular structure (that 
is they have vertices, lines, faces and also 3-cells, 4-cells, and so on up to D-cells), thus they are 
tensor models. 

Definition 2.3. A (D + \)-colored model is a probability measure dv. 

D D 

where 
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• (j/ : ^ C are D + 1 complex (Grassmann) random fields^; 

• : F{X'') F{X') are L> + 1 covariances; 

• K,K : X*^ X • • • X X^ — )• C are two vertex kernels. 

The closed Feynman graphs of these colored models are termed closed (-D + l)-colored graphs. 
One may give a constructive approach to their description but we prefer to state the end result. 

Definition 2.4. A closed (D + l)-colored graph is a graph Q = {V,£) with vertex set V and 
edge set £ such that: 

• V is bipartite, that is, there is a partition of the vertex set V = V U V , such that for 
any element / S £, then / = {v,v} where v £ V and v €z V. Their cardinalities satisfy 
|V| = 2\V\ = 2\V\. 

• The edge set is partitioned into D + 1 subsets £ = lj£o '^here £^ is the subset of edges 
with color i. 

• It is {D + l)-regular (i.e. all vertices are {D + l)-valent) with all edges incident to a given 
vertex having distinct colors. 

We shall call the elements v £ V {v G V) positive (negative) vertices and draw them with 
the colors clockwise (anti-clockwise) turning. The bipartition induces an orientation on the 
edges, say from v to v. As is to be expected, we shall only be dealing in with connected graphs 
in the sequel. Finally, an important point on notation, we denote by i,j,...,k the following 
complement in the set of colors {0, . . . , D} \ . . . , k}. 

Remark 2.1. A colored model with real fields can be defined as: 

D 

Its graphs are colored, but their vertex set is not bipartite. 



3 A brief review of (colored) matrix models 

There are many excellent reviews [53] of random matrix models in the literature. This section is 
not intended to be one. Rather we want to briefly recall the most prominent features of matrix 
models, which will then be generalized one by one to colored tensor models. In the case of 
matrix models the colors do not play an important role, being merely a decoration. However, 
having in mind the generalization to higher dimensions we shall include them in the discussion 
below. 

The particular matrix model in which we are interested is a perturbed Gaussian measure for 
three N x N non-hermitian colored matrices, with colors i = {0, 1,2} (extensively analyzed in 
the literature [40, 51, 52]): 

m9 , Mi , Mi , 

no' ni ' n2 ' 

where no = (n*^^,n^^), ni = (^n^^,n^^^, n2 = {it?^,v?^^. 

We will detail below the independent identically distributed (i.i.d.) colored matrix model, 
having covariance and vertex kernels: 

^^i,ni ^ n ^niin^i ; -^no.nina = ,n3i ) ^hoMim ^ "7^ 11 ^n^^ ,nii ■ 

j+i i<j i<j 



■^The domain of definition X' of the fields need not be the same. 
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The reader can convince oneself that the above is just a rather cumbersome way to describe the 
partition function: 

Z = e^= f Y[[dM\dM^)^]e-^^^'^''''^''''^'^-VM'''^^^^ (3.I) 
which can also be written, after the change of variables M = ^/NT, in the more familiar form: 

^ /"-TTr -si, -Af( E, Tr[T'(T01']+ATr[TOTiT2]+ATr[{TO)t(Ti)1'{T2)t] ) 

i 

We will denote E the free energy per degree of freedom E = N^'^F. 

Graphs: The graphs of a matrix model are ribbon graphs comprising of ribbon lines and 
ribbon vertices. The sides of the ribbons are called strands, and closed strands are called faces. 
In the colored matrix model above, the ribbon lines have a color: 0, 1 or 2. The strands (which 
correspond to the matrix indices n*-') have two colors: ij. For instance, the strand colored by 12 
is shared by the ribbon lines of color 1 and 2 at a vertex. The colors of the strands are conserved 
along the ribbon lines, hence the faces are indexed by couples of colors, see Fig. 1. Due to the 
presence of the colors the ribbon graph representation is, as a matter of fact, redundant. Indeed, 
one can represent a graph using only point vertices and colored lines. Every such graph can be 
uniquely mapped onto a ribbon graph by adding faces for the cycles made of lines of only two 
colors. 

For every graph one can construct a dual triangulation by drawing a triangle dual to every 
ribbon vertex and an edge dual to every ribbon line. The faces of the ribbon graph correspond 
to the vertices of the triangulation, see again Fig. 1. Thus, every graph is dual to an orientable 
surface. The free energy F is the partition function of connected surfaces. 




Figure 1. Colored ribbon graphs. 



Observables: A convenient set of observables in matrix models are the so called loop observ- 
ables. These are traces of products of matrices corresponding to external faces of connected 
open ribbon graphs. The colors impose some restrictions on these observables. As an exam- 
ple, the reader can convince oneself that a matrix M* can only be followed by either (M*)^ or 
by m(*+i hence: 



'TV[MO(MO)t], 
Valid: <^ Tt:[M°M^M'^], 

^Tr[MO(MO)t(Mi)tMi], 



Invalid: Ti[M° {M^)^My M^] 



Amplitudes and the 1/N expansion: As all lines connect a A- vertex with a A-vertex, a closed 
connected graph has an even number of vertices |V| = 2p. The vertices are trivalent, thus 
the number of lines \£\ computes to \£\ = ||V| = 3p. The amplitude of a ribbon graph may 
be readily computed in terms of the numbers |V| = 2p, \£\ and \F\ (the number of faces of 



Colored Tensor Models - a Review 



13 



a graph). Starting from (3.1) (of course, the same result is obtained, starting from (3.2)), we 
obtain a contribution N~^^'^ per vertex and a contribution N per face, hence: 

A{g) = {X~X)PN-P+^ = {X\yN^-^+^ = (AA)PiV2-29, 

with g the genus of the graph. As the genus is an nonnegative integer, the free energy F of 
a matrix model supports a 1/A^ expansion. The leading order graphs are those with genus g = 0, 
that is, planar graphs corresponding to spherical surfaces. Higher genus surfaces are suppressed 
in powers of 1/iV [42, 46, 93]. 

Continuum limit: The planar sector of the 3-colored matrix model can be analytically sol- 
ved [40, 51, 52]. After defining the constant t hy t = (XX)^^^, one can rewrite the derivative of 
the free energy of the model (after a change of variables A = A^/^T) in parametric form as: 

dE 

t^ = U = ^{l-^z), t = z{l-2z), 

which is solved by: 

1 - VT^^ 3 3 1 (1 - 8t)3/2 

z = If = 1 1- 

4 ' 4 ^ 8t 32*2 ^ 32^2 ' 

When t approaches the critical value tc = |) the free energy of the matrix model (i.e. the 
partition function of connected surfaces) exhibits a critical behavior: 

-Ssing ~ {tc — t)^/"^. 

The average area of the connected surfaces, {2p) ~ tdt In £' ~ diverges, hence, large surfaces 
(triangulated by many triangles) dominate. 

Loop equations: The integral on two colors, say 1 and 2 is Gaussian and can be explicitly 
performed to express the partition function of the colored three-matrix model as an integral over 
only one matrix: 

Z = j [dTdT^]e'"'A^'^^'<^-'^^^'^^] = I [drdrt]e-^[^'^]+^-^ ^^^[(^'^)''] . 

The model becomes the most general model for one (non-hermitian) matrix by replacing the 

coefficients by an independent tp for each of the operators Tr[(r^r)^] in the effective 

action for the last color: 

/oo 
^-NTrViT^T)^ y^J^tpiT^Tf. 
p=l 

The Schwinger-Dyson equations (SDE) of this one-matrix model are: 

= J [dM] -J^{[(MM^rMU e-^Tr[v{MMt)]^ 

n 

Y,mM^)%a[{M^Mr-%) - Af( j;jt,-[(MMt)-M],,[Mt(MMt)i-i],„ 

k=0 3=1 
which, summing over a and 6, become: 

n 

l^^i:v[{MM^f]l^[{MM^ff-^^ - N^j tj(l^[{MM^Y+^]^ =0. 

fc=0 3 
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Every insertion of an operator Tv[{MM^y] in the correlation function can be re-expressed as 
a derivative of V{MM^) with respect to tj. Consequently, the SDEs can be written alterna- 
tively as: 

LnZ = 0, for n > 0, 

where 

2d 1 ^ . d 



2d 1 5^ °° 



k=l j=l 

where the derivatives w.r.t. tj, with j < are understood to be omitted. A direct computation 
(involving some relabeling of discrete sums) shows [71] that the -L„'s respect the commutation 
relations of (the positive operators of) the Virasoro algebra: 

[Lm, Ln] = (m - n)Lm+n for m, n > 0. 

Note that as we only deal with Lm, m > 0, we certainly do not obtain the central charge term. 



4 Topology of colored graphs 

Edge-colored graphs have been extensively used in topology [64, 104] to study manifold crystal- 
lization. We shall present in this section the topology of {D + l)-colored graphs, but without 
restricting to manifold topologies. The colors encode enough topological information to con- 
struct a D-dimensional cellular complex, rather than just the naive 1-complex of a graph. An 
intuitive picture of the complex associated to a graph is presented in Remark 4.1. Essential to 
the construction of this D-dimensional graph complex are the d-cells (or d-bubbles as we term 
them), for all d = 0, . . . , D. This nested cellular structure is precisely the information encoded 
in the colors. With this at our disposal, one can subsequently define a rather rich cellular 
(co-)homology upon the graph complex. 

To aid the reader's navigation through this rather technical section, let us give a flavor of 
the subsequent subsections. 

Cellular structure and pseudo-manifolds: We shall define the d-cells of the graphs, which are 
encoded by the colors, and explain their nested structure. The aim of this subsection is to show 
that a {D + l)-colored graph is dual to a D-dimensional simplicial pseudo-manifold. 

Homology and homotopy: The nested structure of the d-cells permits a definition of colored 
boundary operators and their related colored homology groups. We then supply a presenta- 
tion of the first homotopy group. This subsection stands somewhat alone with respect to the 
development of the rest. 

Boundary graphs: Here, we examine open rather than closed {D + l)-colored graphs. These 
open graphs come equipped with a boundary. We show that these boundaries are themselves 
colored graphs, but this time the colors are associated to the vertices. We shall need these 
definitions when we later examine graph factorizations (below) and the embedded matrix models 
(Section 6). 

Combinatorial moves: This and the following two subsections contribute indispensable know- 
ledge for a thorough analysis of the 1/A^-expansion of tensor models (Section 5), their critical 
behavior (Section 7) and the underlying quantum equations of motion (Section 8). We explain 
A;-dipole creation and its inverse (/c-dipole contraction). These are graph manipulations which 
preserve certain combinatorial features of the graph, although not necessarily its topology. This 
allows us to define the notion of combinatorial equivalence. We conclude by highlighting 1-dipole 
contraction, which plays a significant role in Section 5. 
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Jackets and degree: Jackets are Riemann surfaces embedded in a specific way inside the 
(D + l)-colored graphs. We provide the rules for their identification and show that the number of 
such surfaces, for a given graph, is fixed by its dimension D. Just as the scaling of matrix model 
graphs is controlled by the genus of the Riemann surface associated to that graph, we discover 
later that a generalized concept, known as degree, controls the scaling of higher-dimensional 
tensor model graphs. We define the degree in this subsection; it is simply the sum of the genera 
of all the jackets of a graph. 

Topological equivalence: Although not all /c-dipole creations/contractions preserve topology, 
there is a subset that possess this property, that is, they are homeomorphisms. With this in 
hand, we are finally at the stage at which we can provide one of the most important results used 
in the analysis of tensor models, namely, if the degree of a graph is zero, then it is a sphere. 

Graph factorization: We conclude this section by developing partitions on the set of jackets for 
any graph. An element of the partition is known as a PJ-factorization and each one has a pleasant 
property: in order to know the degree of the graph, one needs to know merely the genera of the 
jackets within any one PJ-factorization. Interestingly, this factorization places a spotlight on 
a constraint satisfied by the degree (it cannot take just any value). An understanding of this 
topics is instrumental for grasping the manner in which matrix models are embedded within the 
tensor model structure (Section 6) and for the analysis of classical solutions of tensor models 
(Section 9). 

4.1 Cellular structure and pseudo-manifolds 

At the outset, we need to construct d-bubbles [82, 84] for the graph complex Q. 

Definition 4.1. The d-bubbles of a graph are the maximally connected subgraphs comprising 
of edges with d fixed colors. 

The d-bubbles are denoted by ^S^^^ '*'', where the color indices are ordered ii < i2 < ■ ■ ■ < id 

and p labels the various connected components with the same colors. We denote B^'^ the number 
of d-bubbles of the graph. Note that the 0-bubbles are the vertices of Q, the 1-bubbles are the 
edges of Q. The 2-bubbles are the faces of Q. 



3 3 3 3 

Figure 2. 3-bubbles of a graph in _D = 3. 

A graph in D = 3 and its 3-bubbles are presented in Fig. 2. The 3-bubbles are indexed by 
the colors of their lines, namely from left to right 123, 023, 013 and 012. The 2-bubbles are the 
subgraphs consisting each of two lines, of colors 01, 02, 03, 12, 13 and 23. The 1-bubbles, as 
already mentioned are the lines 0, 1, 2 and 3, while the 0-bubbles are the vertices. 

Now, while the definition of d-bubbles seems sensible, we should check they possess a D- 
dimensional cellular complex structure. Constructing the dual finite abstract simplicial comp- 
lex [100] facilitates this analysis. It goes without saying that the graph complex and the dual 
complex are the same topological space. As a quick note, we write ^ C ^, if ?^ is a subgraph 
of Q. Then, to build up the dual complex [84], we first assemble all the D-bubbles of Q into a 
set A: 




^ = {4p)|Vp, zG{0,...,D}} 
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We consider the subsets of A indexed by the various (D + 1 — d)-bubbles oi Q , d £ {1, . . . , D + 1}: 

Note that for a given bubble ^^^^ each k G {1, . . . , d}, there exists a unique D-bubble 

D ^«i«2---«d^ namely, the maximal connected component (in Q) obtained by starting from 

"^(k)^ and adding lines of all colors except ik- Thus, the cardinality of a is d. Any 

subset T G a is indexed by a choice of subset 5 C {1, . . . , d}: 

and so r = fj , where B^^^ "' is the unique subgraph obtained by adding the 

colors is, s £ S, to the subgraph ^*^*2--.«d^ 

Crucially, the sets a are the {d — l)-simplices of a finite abstract simplicial complex, 

as the collection (multi set) 

A = [a^i^) I eg}, 

is such that VcjGA, rCcT=>rGA. The cardinality of cr S A is d (it corresponds to 
a {D + 1 — (i)-bubble) and so its dimension is d — 1. In fact, since A is non-branching, strongly 
connected and pure, it is a D -dimensional simplicial pseudo-manifold [84].^ 

Remark 4.1. The vertices of the graph correspond to the D-simplices of the simplicial complex. 
The half-lines of a vertex represent the {D— l)-simplices bounding a D-simplex and have a color. 
Any lower-dimensional sub-simplex is colored by the colors of the D — 1 simplices sharing it. In 
Fig. 3 we sketched the dual complex in D = 3 dimensions. The vertices are dual to tetrahedra. 
A triangle (say 3) is dual to a line (of color 3) and separates two tetrahedra. An edge (say 
common to the triangles 2 and 3) is dual to a face (2-bubble of colors 2 and 3). A vertex (say 
common to the triangles 0, 2 and 3) is dual to a 3-bubble (of colors 0, 2 and 3). 




Figure 3. The dual complex in _D = 3. 



^For the sake of self-containment, we provide here a concise explanation of the above comment. A D- 
dimensional svmplictal pseudo-manifold is a finite abstract simplicial complex with the following properties: 

• non-branchmg: each (D — l)-simplex is a face of precisely two D-simplices; 

• strongly connected: any two D-simplices can be joined by a "strong chain" of D-simplices in which each 
pair of neighboring simplices have a common {D — l)-simplex; 

• pure (that is, dimensional homogeneity): each simplex is a face of some D-simplex. 
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Remark 4.2. One can give a representation of a colored graph which includes the strands. 
Recall that the strands represent the faces of the graph and thus are identified by pairs of 
colors. Every line of color i is represented by D parallel strands with colors ij, j ^ i. Inside 
every colored vertex the strand ij is the strand common to the half-lines i and j incident at 
the vertex. Therefore, the ribbon representation of matrix models graphs becomes a stranded 
representation for graphs made of lines and vertices. This is depicted in Fig. 4. 




Figure 4. Stranded representation. 

This representation is often used in the literature. Being redundant and very cumbersome, 
we shall not use it at all in this review. 

Remark 4.3. The graphs of the real colored model are also simplicial pseudo-manifolds. Indeed, 
the fact that the vertex set is bipartite not played an important role so far. This will be the 
case, up to trivial generalizations, for many of the results we present in the sequel. The bipartite 
condition is equivalent to the orientability of the pseudo-manifold [44]. 



4.2 Homology and homotopy 

This section details the topology of the graph complex. It can be skipped at a first read, 
especially by the reader interested in the critical behavior of tensor models, as it does not play 
a prominent role in the latter. The homology of the topological space defined by a graph is 
studied by means of a colored homology defined for the graph complex [82]. 

Definition 4.2. The d'th chain group is the group finitely generated by the d-bubbles: 

2^ %) '^{p) ' %) 

The chain groups define homology groups via a boundary operator, 
Definition 4.3. The d'th boundary operator dd acting on a d-bubble ^jp-)"*'' is: 

• for d > 2, 

q-l n n-.id 

(k) (p) 

which associates to a d- bubble the alternating sum of all {d— l)-bubbles formed by subsets 
of its vertices. 

• for d = 1, since the edges connect a positive vertex v to a negative one v: 

ddBl) =v-v. 

• for d = 0, dov = dov = 0. 
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The colored boundary operators define a homology as d^-i 084 = [82] and thus we define 
the d'th colored homology group to be Ha = ker(9d)/Im(9d+i)- We list some properties of the 
maximal and minimal colored homology groups in the lemma below [82]. 

Lemma 4.1. Denoting the total number of d-bubbles by B^'^ {hence B^^^ its the number of 
vertices, B^^^ is the number of lines, etc.): 

ker(5o) = 0Z, Im(5i) = Z, hence Ho = Z, 

ker(9D) = ^, Ini(9D+i) = 0, hence Ho = TL. 

At the same time, some more useful information falls into our lap: 
Im(5D)= Z, ker(ai)= Z. 

Remark 4.4. A finite presentation of the fundamental group of the graph Q is obtained by 
associating a generator hi to all edges I G G (apart from those edges lying on a tree in Q, which 
we fix to the identity e) and a relation to all faces B^^^^y. 

^B-, = n = 1' 
'^^(^) 

where the product follows the boundary of the face BV-^ in question and e(^) is +1 if our direction 
around the boundary agrees with the orientation of / and —1 if not. One often comes across 
tensor models that weight graphs according to this characteristic. 

Example 4.1. We present some examples of graphs and their associated homology groups. 
The black (white) vertices are the positive clockwise turning (negative anti-clockwise turning) 
vertices. 



,-...o 




Figure 5. Examples of 3-dimensional graphs. 

• Consider the leftmost 4-colored graph in Fig. 5. Its homology groups are: 

Hq = Z, Hi = Z2, H2 = 0, H^ = Z, 

matching those of MP^. In fact, one can show that the graph represents MP^, but this 
requires more effort than merely computing homology groups. 

• A second example is given by the central graph in Fig. 5. Its homology groups are: 

Hq = Hi = Z, Hi = H2 = 0, 
which are those of the 3-sphere §^ (we shall see later that this graph represents a sphere). 
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• For the final example, we take the rightmost graph in Fig. 5. Its homology groups are: 

Hq = 'L, Hi=ZeZ, H2 = 0, = z. 

This graph represents a pseudo-manifold with two isolated singularities. In 3-dimensions, 
the topological singularities arise just at the vertices of the simplicial complex and we find 
in this case that the links^ of two vertices are homeomorphic to the torus. 

4.3 Boundary graphs 

Up to now we dealt with vacuum graphs, that is, closed {D + l)-colored graphs. Non-vacuum 
graphs arise when one evaluates observables and have a natural interpretation of Z)-dimensional 
pseudo-manifolds with boundary. As such they play an important role in colored tensor models. 
The boundary is itself a {D — l)-dimensional pseudo-manifold. To begin, let us describe these 
non- vacuum graphs: 

Definition 4.4. An open {D + \)-colored graph is a graph Q satisfying some additional con- 
straints: 

• It is bipartite, that is, there is a partition of the vertex set V = V UV , such that for any 
element I E £, then / = {v, v} where v £ V and v £ V. 

• The positive vertices are of two types V = Vint + Vd, where T^nt is the set of {D + 1)- 
valent internal vertices and the elements of Vg are 1-valent boundary vertices. A similar 
distinction holds for negative vertices. 

• The edge set is partitioned into D + 1 subsets S = Uiez^+i where f * is the subset of 
edges with color i. Furthermore, each £^ = £^^^U£l^^, such that internal edges £^^^ join two 
internal vertices, while external edges fg^t join an internal vertex to a boundary vertex. 

• All edges incident to a D + 1 valent vertex have distinct colors. 

We consider only connected open graphs, from which we construct the boundary graph [87] 
as follows. 

Definition 4.5. The boundary graph Qq of an open [D + l)-colored graph Q comprises of: 

• the vertex set Vq = VqUVq. We stress that it is not bipartite with respect to this splitting. 
The vertices inherit the color from the external edges of Q upon which they lie, so that 
a more appropriate partition is Vg = UieZo+i where Vq denotes the set of boundary 
vertices with color i. 

• the edge set £q = \Ji:^j^Xn+i^d ^ where = {v,w} G £q exists if there is a bi-colored 
path from v to w m Q consisting of colors i and j. Thus, the lines £q inherit the colors of 
the path in Q. 

These boundary graphs, possess a number of additional properties. The line l^^ £ £q can 
only exist if P^ = {v^,w^} or {v^,w^} or {v^ ,w'''} or {v^,w^}, where v^,w^ £ Vq and £ Vq. 

Each boundary vertex is D-valent and for £ Vq, the incident boundary edges are where 
j ^ i- Several examples for D = 3 are presented in Fig. 6. 

Remark 4.5. As it has colored vertices and bi-colored edges, the boundary graph Qq is a priori 
very different from the initial graph Q. The boundary graph can have several connected com- 
ponents. Each component has a cellular complex structure. For d > 1 the boundary d-bubbles 



The link of a vertex in a simplicial complex is the simplicial complex IkA(l^) — {a E A\V ^ a and VUa € A}. 
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Figure 6. Open 4-color graphs and their corresponding boundary graphs. 



{BgY^^^'^''''^^ are the maximally connected components of Qq formed by boundary vertices f and 

boundary edges where ia,ib,ic £ {h, ■ ■ ■ ,id+i}- Following step-by-step the constructions 
for the vacuum graphs, but taking into account that the boundary d-bubbles have d+1 colors, 
one can show that each connected component of Gg is dual to a simplicial complex (and is 
a pseudo-manifold). In fact, the simplicial complex dual to Gq is the boundary of the simplicial 
complex dual to Q. Consequently, one can study its colored homology. 

A particular subclass of boundary graphs are those obtained from open (D + l)-colored 
graphs, such that all external but none of the internal lines have color D. In this case, the open 
graph possesses exactly one internal D-bubble. Additionally, all vertices in the Qg have the same 
color, D and all edges in Qg correspond to the existence of a bi-colored path in Q joining the 
two vertices. Hence, all edges are of the form l^^ for an internal edge of color k. In fact, if one 
deletes the label D from all the vertices and edges of the boundary graph Qg, one finds that it 
is identical to the graph obtained from Q by deleting the external edges. 

4.4 Combinatorial moves: fc-dipole reductions 

The colored graphs support a class of combinatorial moves, termed /c-dipole moves [85, 86, 90], 
which have a well controlled effect on their bubble structure. In a subsequent section, we shall see 
that under a further assumption, the moves implement homeomorphisms of the topological space 
associated to the graph. These moves are crucial to understanding the structure of arbitrary 
terms in the large N expansion of colored tensor models. However, the leading order of the 
expansion can be well understood without going into the details of this section. 

Definition 4.6. A k-dipole dt is a subset of Q comprising of two vertices v, v such that: 

• V and V share k edges colored by ii, . . . , G Z/j+i; 

• V and V lie in distinct [D + 1 — A;)-bubbles: iJ^^y*''' ^ S^^'^"**. 

We say that dk separates the bubbles -B^^y** and S^^y**. Yet more important is how we 
manipulate the graph structure with respect to these subsets. 

Definition 4.7. The process of k-dipole contraction: 
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• deletes the vertices vi and V2; 

• deletes the edges ii, . . .ik] 

• connects the remaining edges respecting coloring, see Fig. 7. 



D 









k 



k - 1 




Figure 7. The process of /c-dipole contraction/creation in the case ip = p — 1 for all p e {1, . . . , fc}. 

The inverse of /c-dipole contraction is called k-dipole creatioriJ . As a point of clarification, we 
note that it is certainly not the fc-dipole creation results case that selecting arbitrarily D -\-l — k 
distinctly colored strands, snipping them and inserting the partial subgraph on the left hand side 
of Fig. 7. It is rather the case that we much carefully choose the strands such that the resulting 
{D + \ — A;)-bubbles are distinct. We denote Q /dk the graph obtained from Q by contracting dk- 

Definition 4.8. Two graphs as said to be combinatorially equivalent, denoted if they are 
related by a sequence of /c-dipole contractions and creations. 

In particular, Q Q /dk- 

We now analyze the effect a A;-dipole contraction has on the graph. Consider a vacuum 
graph Q and denote by Vg, £g, Tg its vertex, line and face sets respectively, such that |Vg| = 2p. 
Under a fc-dipole contraction, Q — t- Q jd^ and the cardinalities of their various sets are related 



In words, the number of vertices decreases by two and the number of lines by Z? + 1 (one per 
color). For the faces, note that prior to contraction, they come in four types: the k{k — l)/2 
faces that are comprised exclusively of lines in the dipole; the k{D + 1 — A;) faces that contain 
exactly one line in the dipole; the (Z) + 1 — A;)(Z) — A;) faces that do not contain any line in 
the dipole, but contain one of the two vertices v or v; and the rest, which do not concern us 
here. After contraction of the A:-dipole, the faces made up exclusively by lines in the dipole are 
erased, while the faces containing either v ox v but no line in the dipole are merged pairwise. 
The k{I) + 1 — fc) faces containing exactly one line in the dipole will lose one line, but their 
number is conserved. 

In subsequent sections, 1-dipole contraction will play a significant role and so let us examine 
this case in more detail. Consider a 1-dipole d\ of color i. We can see from above formulae that 
Qjdx has two less vertices. Interestingly, Qjdx also has one less D-bubble than Q of colors i, 
while the number of -D-bubbles of colors j, j 7^ i is left unchanged. Thus, the quantity p — iS'^^ 
is conserved by 1-dipole contraction. Also, the connectivity of the graph remains the same. So 
now let us contract a maximal number of 1-dipoles in Q so that it reduces to some graph Qj with 
2pf >2 vertices. With a little inspection, one notices that Qf possesses exactly one D-bubble 
for each color i, and so, Bq''^ = D + 1. We have arrived at the important inequality: 



by: 



\Vg/d,\ = \Vg\-2 



^ \£g/d,\ = \£G\-{D + l), 
k{k-l) {D + l-k){D-k) 
2 2 



(4.1) 



p- 




Pf 





(4.2) 



'^The contraction/creation of a _D-dipole is superfluous and can always be traded for the contraction/creation 
of a 1-dipole. 
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4.5 Jackets and degree 

In order to gain a better understanding of colored graphs, one would like to define some simpler 
graphs which capture only some of the information encoded in the colors. A first class of sub- 
graphs is already available: the bubbles. The bubbles are, however, themselves colored graphs. 
Hence although they have less colors, they are still relatively difficult to handle. A second class 
of simpler graphs is given by the jackets [23, 85, 86, 90]. The main advantage of the jackets is 
that they are merely ribbon graphs (i.e. they comprise just of vertices, edges and faces), like the 
ones generated by matrix models graphs. The jackets contain all the vertices and all the lines 
of Q but only some of its faces. As they are needed in order to define the 'degree', they play 
a crucial role in the large N expansion of colored tensor models. As we shall see, the jackets are 
Riemann surfaces embedded in the cellular complex. 

Definition 4.9. A colored jacket ^7 is a 2-subcomplex of Q, labeled by a (Z) + l)-cycle r, such 
that: 

• J' and Q have identical vertex sets, Vj = Vg; 

• J and Q have identical edge sets, £j = £g; 

• the face set of J" is a subset of the face set of Q: Tj = {/ G Tg \ f = (r''(0), r'?+^(0)), q G 
Zd+i}. 

It is evident that and G have the same connectivity. In actual fact, a given jacket is 
independent of the overall orientation of the cycle, meaning that the number of jackets is in 
one-to-two correspondence with (D + l)-cycles. Therefore, the number of independent jackets 
is D\/2 and the number of jackets containing a given face is {D — 1)!.^ 

The jacket has the structure of a ribbon graph. Note that each edge of lies on the boundary 
of two of its faces, thus corresponds to a ribbon line in the ribbon graph. As we said, the ribbon 
lines separate two faces, (r~^(i),i) and (z,r(i)) and inherit the color i of the line in J. Ribbon 
graphs are well-known to correspond to Riemann surfaces, and so the same holds for jackets. 
Given this, we can define the Euler characteristic of the jacket as: x{J^) = \ ~ \^j\ + = 
2 — 2gj^ where gj is the genus of the jacket^. 

In D = 2, the (unique) jacket of a (2 + l)-colored graph is the graph itself. In D = 3, an 
example of a graph and its jackets (and their associated cycles) is given in Fig. 8. For instance 
the leftmost jacket corresponding to the cycle r = (0123) contains only the faces 01, 12, 23 
and 30. 




T = (0123) T = (0213) r = (0132) 



Figure 8. A vacuum 4-colorcd graph and its jackets. 

For a (D + l)-colored graph its -D-bubbles are D-colored graphs ^"^i^py Thus, they also 
possess jackets, which we denote by J]^py It is rather elementary to construct the Jl^^-^ from 

®It is, however, sometimes more transparent to over count the distinct jackets by a factor of two associating 
them one to one with cycles. For example, on can count that from the D\ cycles of D + 1 colors, (f — 1)! will 
contain the pair ij and (D — 1)! the pair ji. 

momentary reflection reveals that the jackets necessarily represent orientable surfaces. 
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the J . Let us construct the ribbon graph J'' consisting of vertex, edge and face sets: 

Vj- = Vj, £j- = £j \ £\ Tj- = [Tj \ {(T-i(i),.), (i,r(z))}) U {(r-i(i), t(0)}, 

that is having all vertices of all lines of Q of colors different from i and some faces. Given 
that the face set of J is specified by a (D + l)-cycle r, the first thing to notice is that the 
face set of J'^ is specified by a L'-cycle obtained from r by deleting the color i. The ribbon 
subgraph J'^ is the union of several connected components, Jl^py Each J^^^ is a jacket of the 

Z?-bubble -B^p)- Conversely, every jacket of is obtained from exactly D jackets of Q. To 

realize this, consider a jacket Jl^py It is specified by a D-cycle (missing the color i). On can 
insert the color i anywhere along the cycle and thus get D independent [D + l)-cycles. 

More generally, the d-bubbles are d-colored graphs and they also possess jackets which can 
be obtained from the jackets of Q. 

Consider once again Fig. 8. Applying our procedure to the jacket (0123) leads to the three 
jackets (123), (023) and (012). Each of these jackets corresponds to a bubble of Fig. 2 and is 
a 2 + 1 colored graph. 

Definition 4.10. We define: 

• The convergence degree (or simply degree) of a graph Q is a;(^) = ^jQj, where the sum 
runs over all the jackets J oi Q. 

• The degree of a fc-dipole dk separating the bubbles ^^J^y** and denoted uj{dk), is 
the smallest of the two degrees ui&f^^^^'') and a;(;B^Jjy*'°). 

The degree is a number one can readily compute starting from a graph. As it will play a major 
role in tensor models, we list below a number of properties of the degree. For understanding 
the leading order in the large expansion one needs to recall from the reminder of this section 
that the degree is a positive number, and equation (4.3) relating the number of faces of a graphs 
with its degree. 

First, the degree of a graph and of its bubbles are not independent. 
Lemma 4.2. The degrees of a graph ijo{Q) and its D-huhhles u}{B^^^-^) respect: 



\p)>- 

i;p 

In particular, by (4.2), p+D-B^^^ > 0, thusuj{g) = ^ ^i^{p)) = Vi,p andp+D-B^^^ = 0. 

Proof. Consider a jacket J dQ. The number of vertices and lines of J are: = |Vg| = 2p 
and iVj-j = \Eg\ = (Z) + l)p, respectively. Hence, the number of faces J is = {D — l)p + 
2 — 2gj. Taking into account that G has jackets and each face belongs to {D — 1)! jackets: 

, , D(D-l) 2 v-^ D(D-l) 2 , , 

Each of the D-bubbles (with \ = ^P\p^) is a D-colored graph and thus an analogous 
formula to (4.3) holds for each D-bubble: 
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Each vertex of G contributes to D + 1 of its D-bubbles and each face to D — 1 of them. Thus, 
J2rpP\p) = {D + l)p and Xli p l-^gi I = (-^ ~ 1)1-^^1- Summing over D-bubbles in (4.4) and 
dividing by D — 1 yields: 

One equates (4.5) to (4.3) to prove the lemma. ■ 

Second, due to (4.1), the degree of a graph changes under a fc-dipole contraction. 
Lemma 4.3. The degree of Q and Q/dk are related by: 

= ^^^-^{{D + l)k -k^-D)+ u:{g/dk). 

One can track in more detail the effect of a 1-dipole contraction on the degrees of the D- 
bubbles of a graph. Let us denote by di a 1-dipole of color i separating the two D-bubbles 

and B'^(^py For j / i or {j = i, p a, (3} we denote B-j^^^ (resp. B-'^^-^/di) the D-bubbles before 

(resp. after) contraction of di. The two Z^-bubbles B^^^ and are merged into a single bubble 

'^(a)u(/3)/^i after contraction. 

Lemma 4.4. The degrees of the bubbles before and after contraction of a 1-dipole respect: 

'^(^(p)M) = ^(^(p))' {« = i, P/a,/3}, 

^(%)u{«M)='-(4))+^(%))- 

Proof. Only the D-bubbles containing one (or both) of the end vertices v and v of di are 
affected by the contraction. 

For j 7^ i, using formulae analogous to (4.1) for D-colored graphs, one sees that any jacket 

•J^p^/di of B'l^^^/di has 2 vertices less, D lines less and D — 2 faces less than the corresponding 
jacket ^/ N of Therefore, 07, =91 and the degree is conserved. 

Any jacket <^(q)u(/3)/'^i of the -D-bubble 'S(q)u(/3)/^1' obtained by gluing two jackets J7^*^^ and 
J7(*^j of and -S^^^. It follows •^(q)u(/3)/'^i ^ vertices less, D lines less, D faces less and 1 



connected component less than the two jackets and J'L-^- Hence g ~ = 9 ^1 + 9 , 



and the lemma follows. 

An important consequence of Lemma 4.2 is that two graphs related by a 1-dipole contraction 
have the same degree. 

4.6 Topological equivalence 

We now include the topology in the picture. We shall utilize a fundamental result from combi- 
natorial topology [64, 104]: 



Theorem 4.1. Two pseudo-manifolds dual to Q and Q/dk are homeomorphic if one of the 
bubbles ^*°y*''-i or B'^^^j'^'"'^ separated by the dipole is dual to a sphere S^~^. 

This allows us to propose another equivalence relation on the set of colored graphs. 
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Definition 4.11. Two graphs, Gi and Gf, are said to be topologically equivalent, denoted if 
they are related by a sequence of dipole contraction and creation moves satisfying an additional 
property: for any dipole move in the sequence, at least one of the bubbles separated by the 
dipole is a sphere. 

It is in principle very difficult to check whether a graph is a sphere or not. However we 
can establish the following partial result, crucial for the large expansion, including at leading 
order, of colored tensor models. 

Lemma 4.5. If uj{G) = then G is dual to a sphere . The reciprocal holds in D = 2. 

Proof. We use induction on D. In Z) = 2, ^ is a ribbon graph and its degree equals its genus. 
For D > 3, as u}{G) = 0, Lemma 4.2 implies that uj(B^^^-^) = and, using the inductive hypothesis, 

all the bubbles ^B^^^ are dual to spheres S^~^. Any 1-dipole di will separate a sphere, so that 

G/di G and by Lemma 4.3, uj{G/di) = u}{G) = 0. We iteratively contract a full set of 
1-dipoles to reduce ^ to a final graph Gf, such that Gf G, ^{Gf) = and Gf does not posses 
any 1-dipoles. It follows that Gf has B^^^ = D -\- 1 remaining D-bubbles (one for each colors i) 
and, by Lemma 4.2, pf + D — B^^^ = 0. We conclude that pf = 1 and Gf represents the coherent 
identification of two D-simplices along their boundary i.e. it is dual to a sphere . ■ 

A last result we shall need is a lower bound for the degree of a graph as a function of the 
degrees of its -D-bubbles with fixed colors . 

Lemma 4.6. Let G be a D + 1 colored graph and B^^^ its D-bubbles with colors D. Then 

p 

Proof. Consider a jacket J of G- By eliminating the color D in its associated cycle we obtain 
a cycle over 0, . . . , — 1 associated to a jacket J'^^ for each of its bubbles. As ribbon graphs, jTj-^-j 
are in one-to-one correspondence with disjoint subgraphs of J. One obtains these subgraphs by 
deleting the lines of color D and joining the strands {7r~^{D),D) and {D,-k{D)) in mixed faces 
corresponding to {■k~^{D),7:{D)) in J'^^ [90]. Consequently: 



jD ■ 

P 



As we already mentioned, every jacket J'P. is obtained as subgraph of exactly D distinct 



'(p) 

jackets (corresponding to inserting the color D anywhere in the cycle associated to >J^^)- 
Summing over all jackets of G we obtain: 

J P ^"^ 

4.7 Graph factorization 

Earlier, we introduced jackets, which are instrumental in the definition of the convergence degree. 
Later, we shall see that the amplitudes of many colored tensor models capture information 
about these jackets rather than higher- dimensional subspaces. Let us recall that the tensor 
models generate Feynman graphs using identical building blocks, namely the interaction vertices. 



26 



R. Gurau and J. P. Ryan 



Interestingly, it emerges that there is a seed for the jackets directly within this building block 
and thus in the tensor model action itself. Identifying this basic structure will be our purpose in 
this section; we shall utilize it later to provide an interesting perspective on the colored tensor 
models in question. 

Our first port of call is the interaction vertex. We usually picture this as a single vertex with 
D + 1 distinctly colored half-edges emanating from it. Obviously, this is a viable graph in its 
own right for the interior of a manifold with boundary. It is more convenient here to treat it as 
such and to formulate our analysis in terms of its boundary graph. Using the rules set out in 
Section 4.3, we note that the boundary graph has D + 1 vertices, each of which is connected to 
all the others by a single boundary edge. In other words, the boundary graph is the boundary 
graph of a D-simplex. On top of that, each of the boundary vertices possesses a distinct color 
(inherited from the corresponding colored half edge of the interior), while each boundary edge 
is labeled by two colors (inherited from the corresponding interior face). 

The introduction of this boundary graph allows us to perform a more concise investigation. 
Its power lies in the fact that any result found for the boundary graph can be easily translated 
back to the interior since, in this case, the vertex and edge sets of the boundary are in one-to-one 
correspondence with the edge and face sets of the interior, respectively. 

To keep the following manipulations as succinct as possible, we shall also need some elemen- 
tary definitions from graph theory. 

Definition 4.12. Some graph theory definitions: 

• A complete graph is one in which every pair of distinct vertices is connected by a unique 
edge. Moreover, we shall be interested in complete graphs with distinctly colored vertices. 
We shall denote such a colored complete graph with n vertices by /C„. 

• A k-f actor Ik is a k-regular spanning subgraph oiQ. In other words, its vertex set coincides 
with that of G and all vertices are A:-valent. 

• A k- factorization of Q partitions the edge set of Q into disjoint fc-factors^*^ . 

• The graph subtraction operation is based on the set-theoretic one. In particular, we shall 
be interested in Q — Ti^ the graph Q with the edges of a 1-factor removed from its edge set. 

The boundary of a colored D-simplex is the colored complete graph 1Cd+i- Furthermore, 
a 2-factor corresponds to a (D -|- l)-cycle in the vertex set of the boundary that is a cycle in 
the set of colors. Thus, on construction of the Feynman graph this 2-factor determines a jacket. 
In a moment, we shall also need information about the face subset of a generic tensor model 
graph Q that is determined by a 1-factor of K,d+i- For this, we shall introduce a new object 
in called a patch, corresponding to the image of this 1-factor. 

Definition 4.13. A colored patch "P is a 2-subcomplex of for odd D, labeled by a 1-factor Xi 
in K-D+i-, such that: 

• V and Q have identical vertex sets, V-p = Vg] 

^°For small values of k, these fe-factors often come under different names. 

— A matching in a graph is a set of edges without common vertices. A perfect matching is one which contains 
all vertices of the graph. Thus, a perfect matching and a 1-factor are equivalent concepts. 

— A cycle is an alternating sequence of points and edges, {wo, eoi, ^i, • • • , Wn, e„o}, such that eu+i = ViVi+i 
and all vertices are distinct. We denote it by r = («o, ui, • • • , v„), leaving the edges implicit. A Hamiltonian 
cycle is a cycle which contains every vertex of Q and thus it is a 2-factor. 

— An edge decomposition of a graph is a partition of its edges into subgraphs. A Hamiltoman decomposition 
is an edge decomposition consisting of Hamiltonian cycles and so it is a 2-factorization. 
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• V and Q have identical edge sets, E-p = £g; 

• the face set of "P is a subset of the face set of G: T-p = {f & J^g \ f = (ij) S Ii}. 

Unlike the jackets, a patch does not correspond to a Riemann surface embedded in the cellular 
complex. On the contrary, any two faces in a patch are either disjoint or intersect in a finite 
number of (zero-dimensional) points. 

Before going any further, we recall a rather well-known result that subsequently plays a sig- 
nificant role. It relates to a permissible factorization of /C_d+i and as such, it is the simplest and 
most generic. 



Proposition 4.1. Say n G N. There exists a 2-factorization of fC2n+i dnd fC 



2n+2 



Proof. We employ an explicit construction due to Walecki [1] in order to prove this proposition. 
For /C2n+i; label the vertex set by Z2n+i and construct the 2- factor: 

(2n, 0, 1, 2n - 1, 2, 2n - 2, 3, 2n - 3, . . . , n - 1, n + 1, n). 

Next act on this cycle with the permutation p = (0, 1, 2, . . . , 2n — l)(2n). This generates another 
2-factor disjoint from the first. Repeating this action another n — 2 times generates n disjoint 
2-factors and thus a 2-factorization. Perhaps an illustration for n = 3 shall serve to clarify 
matters, see Fig. 9: We form a 2n-gon with the vertex labeled by 2n at the center. We can think 
of the action of p as rotating the labels (passive) or equivalently rotating the cycle (active). 





4 3 4 3 4 3 

/Ct (6,0,1,5,2,4,3) (6,1,2,0,3,5,4) (6,2,3,1,4,0,5) 
Figure 9. 2-factorization of /Cy. 
For fC2n+2, we label the vertex set by Z2„+2 and construct the 2-factor: 

2n -Fl, 0, 1, 2n - 1, 2, 2n - 2, 3, 2n - 3, . . . , -, 2n, y, . . . , n - 1, n -Fl, n 



for n even, 



2n + 1, 0, 1, 2n - 1, 2, 2n - 2, 3, 2n - 3, . . . , 



3n + 1 „ n + 1 

2n, 



. . . ,n — l,n + l,n 



for n odd. Once again by applying p = (0,1,2, ... ,2n — l)(2n)(2n + 1), one generates n disjoint 
2-factors. But in this case, it does not exhaust the edge set of JC2n+2- There is a 1-factor left 
over: 

Ii = {{2n,2n + l)}U{{x,x + n) : X eZn}. ■ 

Note that in our case, the vertices of /Cd+i are distinguished by a color. Therefore, the 
2-factors of this labeled K^+i are in one-to-one correspondence with the jackets of a (-D -|- 1)- 
colored graph Q, that is, there are ^ of them. For odd D, we shall also need to take into account 
the patches. For a labeled /Cd+i, there are DM 1-factors and thus the same number of patches. 

Let us use the factorization result above to formulate an idea of face-set factorization for 
a graph Q. 
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Definition 4.14. A P, J- factorization VJg of a (D + l)-colored graph Q is: 

• for even a subset of its jackets corresponding to a 2-factorization of fCo+W 

• for odd D, a patch corresponding to a 1-factor Xi, along with the subset of its jackets 
corresponding to a 2-factorization of fCo+i — X\. 

A PJ-factorization has the useful properties that all its elements have disjoint face sets, but 
yet it contains all the faces of Q . Note that for even D, it contains ^ jackets, while for odd D, 
it contains ^^-^ jackets and a patch. 

Another result follows for a labeled D-simplex /C^i+i: 

Lemma 4.7. Consider the complete graph fCo+i with distinguished vertices. One may place 
a partition on the set of 2-factors of 1Cd+i such that: 

• For even D, one may partition the set of all 2-factors into (D — 1)1 2- factorizations. 
Moreover, there are D + 1 such partitions. 

• For odd D, one may partition the set of all 2-factors into D[{D — 2)!] 2-factorizations. 
Once again there are D + 1 such partitions. 

Proof. For even D, consider the 2-factorization constructed in Proposition 4.1. We work from 
the passive stance, where the cycle through the vertices is fixed and the labels are permuted. 
In Proposition 4.1, we started from an initial cycle and considered the action of cyclic subgroup 
of the permutation group Sd+i generated by the element {0,1,2, D — 1){D). This subgroup 
fixes the label D and it has D elements (although due to a symmetry in the initial 2-factor, the 
action of this subgroup generates just D/2 distinct 2-factors). There are {D + 1)[{D — 1)!] left 
cosets with respect to this subgroup, each of which corresponds to a 2-factorization of /C^+i- 
Of these left cosets, {D — 1)! fix a given label, say D, and all contain distinct 2-factors. Thus, 
this is a partition on the set of 2-factors and there are D -\- 1 such partitions. 

For odd D, the situation is bit more involved but one follows essentially the same argument. 
Here, there are two vertex labels fixed and a cyclic subgroup of Sd+i is used to generate the 
rest of the 2-factors. This subgroup has D — 1 elements. Then, the {D -\-l)D[{D — 2)\] left cosets 
are such D[{D — 2)!] of them fix a given label, say D, and contain distinct 2-factors. So, once 
again one has a partition on the set of 2-factors and there D -\- 1 such partitions. 

Note that for odd D, the left-over 1-factor is not the same for each 2-factorization, but it is 
uniquely determined once one has all the 2-factors in the factorization. ■ 

Translating this to a graph Q, we get analogous statements for its set of jackets and patches. 

Proposition 4.2. Consider a {D + l)-colored graph Q . One may place a partition on the set of 
jackets of Q such that: 

• For even D, one may partition the set of all jackets into [D — 1)\ PJ-factorizations. More- 
over, there are D -\- 1 such partitions. 

• For odd D, one may partition the set of all jackets into D\{D — 2)!] PJ-factorizations. 
Once again there are D + 1 such partitions. 

What is of more interest to us here is a definition of the Euler characteristic of a PJ- 
factorization. For even D, the factorization contains only jackets, which are Riemann surfaces 
and thus there is an natural extant definition. For odd D, one must extend this to deal with 
patches. We shall now define the Euler characteristic of a PJ-factorization: 



x{VJg) = { 




\£v\ + I-^pI 



(4.6) 



X{J) + X{V) for odd A where x{V) = \\Vv\ 



\\£v\ + \y'vV 
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Note that xCPJ^g) is an integer and xi'P'Jg) = D — (^^'^iy w{Q). This puts some constraint 
on the values of w{Q): 

w{y) = (an integer). 

Additionally, we notice that to obtain the degree of a graph we need only topological informa- 
tion attached to a subset of the jackets (and perhaps a patch), that is, those in a single PJ- 
factorization. This definition reveals a nice view of Lemma 4.2 relating the degrees of a graph 
and its D-bubbles: 

xiVJg) = -P+^^T.^i'P'^BlJ^ (4.7) 

which stresses the nested structure existing among the various bubbles of a graph. As a word of 
caution, however, we mention that a single PJ-factorization of a graph does not automatically 
provide a PJ-factorization for all its Z?-bubbles. In passing from a graph to one of its Z)-bubbles, 
one deletes a given color. This color gets deleted from the cycles determining the jackets in the 
PJ-factorization. These cycles are no longer a factorization of the ICd generated from ICa+i by 
removing the vertex and edges containing the chosen color. 

Let us focus briefly on the case when D = 3, where we can readily utilize this machinery to 
uncover some further results about the topology of the 4-colored graphs. 

Proposition 4.3. If D = 3 and Q possesses a spherical jacket then Q is spherical. 

Note that one can give several alternative proofs of this result. 

Proof. If w{Q) = 0, we may apply Lemma 4.5 and we are done. Therefore, we consider the 
case where w{Q) / 0. 

Since Q possesses a spherical jacket, there is an ordering of the colors around the vertex that 
provides a planar representation of the graph. Deleting a single color to obtain the 3-bubbles 
of preserves the planarity of this representation. Thus, each 3-bubble possesses a spherical 
jacket. Since 3-bubbles are Riemann surfaces, they are spheres and w{B'^^p^) = 0. 

Since all 1-dipoles separate spheres, we may iteratively contract a full set of 1-dipoles to 
get a graph Qi Q such that Qi has just four 3-bubbles, one for each color. Now, pick the 
PJ-factorization of Qi that contains the spherical jacket. Then, xi'PJ^gi) = 2 + xi'P) ^-^id along 
with the formula (4.7), this gives us: 

1 ^ 

2 + \J^p\-p = 2 + x{V) = -p+-y2x{'PJT^i ) = -p + 4 =^ |Jp| = 2. (4.8) 

Let the spherical jacket be labeled by the cycle (abed). Then, equation (4.8) states that after 
all 1-dipoles are contracted, we are left with a graph having one face of type (ac) and one of 
type (bd). As an illustration of the planar jacket helps at this point, see Fig. 10. 

We have drawn explicitly the face of type (oc) and left most of the lines of color b and d 
implicit. We shall now show that we can iteratively contract a set of 2-dipoles until we arrive 
at the spherical graph with two vertices, denoted Gf. 

Planarity of the jacket implies that of the faces of types (ab) or (be), there are at least two 
which have only two edges. We have drawn one of these, of type (ab), at the top of the figure. 
Showing the lines d emanating from the vertices v and v, one can see that the face {ab) separates 
two distinct faces of type (cd) and so it is a 2-dipole. The only catch would be if u; = v. In that 
case, the graph must he Gf, else it would have more than one face of type (bd). We contract this 
2-dipole and arrive at ^/ or at a graph with the same properties as Gi- Therefore, it contains a 
2-dipole and we may iterate this procedure until we arrive at Gf - To conclude, G is spherical. ■ 
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We present here another result which relates the jackets to important objects in 3-dimensional 
topology: 

Definition 4.15. A Heegaard splitting of a compact connected oriented 3-manifold Ai is an 
ordered triple (T,,'Hi,T-12)m consisting of a compact connected oriented surface S and two 
handlebodies Tii, T-L2 such that dT-ii = d7i2 = S. S is known as the Heegaard surface of the 
splitting. 

Proposition 4.4. If Q is a manifold, then its jackets J are Heegaard surfaces. 

In fact, this result (see [119]) provides an alternative route to proving Proposition 4.3. If Q 
possesses a spherical jacket, then all its 3-bubbles are spheres. This implies that it is a manifold 
(rather than a pseudo- manifold) and its jackets are Heegaard surfaces. A well-known result in 
3-dimensional topology is that a 3-manifold possessing a spherical Heegaard surface must be 
a sphere [91]. 

5 Tensor models 

We have discussed so far the topological and cellular structure of colored graphs. These graphs 
arise due to the coloring of the fields, and thus the results hold for (J) + l)-colored tensor models 
irrespective of their precise details, that is, the fields, their arguments, the vertex kernels, the 
covariances. From now on, we shall start making specific choices for the action of our tensor 
model. The amplitudes of graphs and the physical interpretation of our results depend strongly 
on the particularities of the model. 

Once again, we present a digest of forthcoming subsections. First and foremost, this section 
as a whole details the expansion(s) appropriate for colored tensor models. 

Introducing colored tensor models: We introduce the independent identically distributed 
(i.i.d.) probability measure for colored tensor fields upon which we perform essentially all our 
subsequent analysis. The reader should keep in mind that similar results hold for many other 
choices of measure. We will present in the conclusion a more detailed discussion on the degree 
of generality of our results. 

Boundary graphs as observables: We mention as a quick note how the boundary graphs of 
Section 4.3 arise as moments of this probability measure. 

Amplitude: We write down the amplitude associated to each graph in the i.i.d. model. We 
highlight the importance of choosing the correct scaling (with respect the to parameter A^) for 
the coupling constants of the model. We note that the amplitude of a graph depends solely on 
its degree and moreover, in such a way that the leading order graphs are necessarily Z)-spheres. 

Combinatorial 1/N expansion: In this subsection, we introduce the first of two 1/A^ expan- 
sions. The idea is that for a given number of vertices, one can identify a finite set of graphs, 
combinatorial core graphs, which are 'simplest' with respect to a combinatorial criterion. Every 
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core graph indexes an infinite class of graphs, all having the same amplitude, thus represents 
a term in the 1/A^ expansion of the free energy of our model. Core graphs at higher order (i.e. 
with more vertices) are suppressed in powers of 1/A^. 

Topological 1/N expansion: Here, we repeat the process of the previous subsection with the 
restriction that we are only allowed to use topology-preserving moves. This increases the number 
of core graphs in the topological case and leads to topological core equivalence classes according 
to which the expansion is ordered. 

Redundancies in the 1/N expansion: Despite the many useful properties of this topological 
expansion, it is not as powerful as in the 2-dimensional case, mainly because the i.i.d. amplitude 
is not sensitive to many aspects of higher-dimensional topology. Since it merely depends on the 
degree of the graph, it only captures combinatorial properties. We outline specific superfluities 
of the topological expansion in the i.i.d. scenario. 

First terms of the 1 /N expansion: We present some explicit examples of core graphs at low 
orders. 



5.1 Introducing colored tensor models 

We will always chose the domain of definition of our random functions to be (several copies 
of) a compact Lie group G. Generally, however, it is more convenient to work with Fourier 
transformed fields, that is, fields dependent on the (discrete) representation space of the Lie 
group in question. We will denote by a (large) cutoff in the representations, the range of the 
discrete indices. Thus, generically our colored field is a colored random tensor of some rank r, 
with fii = {n'^i, . . . jU"""}. The crucial feature of the colored models is that one can assign 
weights to the Feynman graphs according to their cellular structure. One might wish, perhaps 
capriciously, to place a weight just on the d-bubbles of colors zi, . . . , i^. This can be achieved 
simply by assigning an index ra*^'"*"^ (resp. ^i^'"*'*) to the field 0*'-' (resp. 0*'-') for k G {1, . . . ,d} 
and utilizing the following covariances and vertex kernel: 



C' = l for i^{ii,...,id}, C\^...,^ = 5 for 



and 



N 



n=l k=l " 

More generally, one can weight several cells (either of the same dimension or not) at the same 
time, use different cutoffs for the various indices and so on. 

We shall concentrate in the sequel on a particular model which assigns equal weights N to 
all the faces (two cells) of our graph. It is the straightforward generalization of a matrix model 
to higher dimensions. Since a line belongs to D faces, our tensors are rank D tensors and we 
denote them (f)^^_, where fij = (n**~^, . . . , n^^ , n*^, . . . , n""'"-'^), with the index n^^ associated to 
the face with colors i and j. The tensors have no symmetry properties in their indices, that 
is, , are independent random variables for each choice of their indices. Furthermore, we 
chose a trivial covariance and vertex kernel: 

^'n, = n fi'' ' . . .no = jyD{D-l)/4: 11 ' 

This is the independent identically distributed (i.i.d.) colored tensor model in D dimensions [85, 
86, 90]. The probability measure we shall deal with in the sequel is therefore: 



2-K 
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D ^ D D 

5(0, ^) = ^Y1 ^\M^nA. + NDiD-l)/4 ^11'^^ + ND(D-1)/A E 11 4 - (^'l) 
«=0 Hi,ni {n} i=0 {n} «=0 

where denotes the sum over ah indices riij from 1 to N. Note that rescahng = N~^/'^Tt 
leads to: 

\i=0 n {n} i=0 {n} i=0 

The scahng of the couphng constant with respect to the large parameter A'^ in (5.1) will be 
explained in Section 5.3. We shall discuss at length the physical interpretation of this model in 
Section 7.4. The notions required for the developing the leading order in the large N expansion 
are the bubbles, the degree, and the relation between the number of faces of a graph and its 
degree in equation (4.3). 



5.2 Boundary graphs as observables 

Tensor models with face weights possess a very convenient property: their observables are in- 
dexed by the boundary graphs. More specifically, the observables are the connected correlations: 

/ 6^ ■ ■ ■ d)U • • • ^ ) . 

As the indices are conserved along the faces, an observable is non-zero if and only if every 
index n'^*^ is paired either with an index n^'^^ or with an index n'^'i. Representing each field 
as a (colored) boundary vertex, and each such identification as a (bi-colored) boundary line 
represents the observable as a boundary graph. Conversely, starting from a boundary graph 
one builds an observable by considering a product of fields, one for each boundary vertex, with 
indices identified according to the boundary lines. 



5.3 Amplitude 

Consider an arbitrary connected D + 1-colored graph Q with amplitude specified by the i.i.d. 
model. Notably, the 6n,n functions associated to the boundary of the faces compose with each 
other until one is left with a exactly one redundant summation per face. Thus, given the scaling 
of the coupling constant in (5.1), the amplitude for a given Q is: 

A{G) = {\~X)P ivl-^el-p^^^ = (AA)P iV^-(iJ^'^(^), 

where we used (4.3) to arrive at the second equality. For the i.i.d. model, we conclude that 
the convergence degree u}{G) indexes the behavior of the amplitude with respect to the large 
parameter A^.^^ In particular, this justifies the scaling of the coupling constant in (5.1). Indeed, 
as we shall see in the sequel, there exists an infinite family of graphs of degree 0. If one does not 
use the appropriate rescaling of the coupling constant, either this family of graphs is increasingly 
divergent at increasingly higher orders (hence the theory makes no sense), or it, along with all 
other graphs, is increasingly suppressed at increasingly higher orders (and in this case the theory 
is trivial, having at most a finite number of divergent graphs. The only scaling which leads to 
a sensible non-trivial theory is the one chosen in (5.1), for which one has infinite families of 
uniformly divergent graphs. 



^'^Hence, the degree acts in higher dimensions as the genus did in matrix models. 
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The perturbative series pertaining to free energy and correlation functions of the i.i.d. model 
are an expansion in the degree: 

F(A,A) = ln(Z) = j;cH(A,A)iV^-(i^", d'^](A,A)= nxfol/^^ 

with s{Q) some symmetry factor. The coefficient C^'^^X, A) is the sum over the set of graphs of 
degree w. The fundamental difference between the two-dimensional case and the general case 
is that, whereas the genus of a surface is a topological invariant, the degree is not in higher 
dimensions. This leads to two distinct expansions, the combinatorial 1/N expansion and the 
topological 1/N expansion, both of which we shall present in detail in a moment. Crucially, the 
leading order in both expansions is given by graphs of degree 0, which we know, from Lemma 4.5, 
are spheres. 

The reminder of this section is dedicated to establishing the full series indexed by the degree. 
For the purpose of the critical behavior in the large N limit however this study is not immediately 
important, as in the large limit all but the graphs of degree zero are suppressed. 

5.4 Combinatorial expansion 

The aim is to give an alternative characterization of Ct'^](A, A) that can be used to index the 
series. From Lemma 4.3, the degree of a graph is invariant under arbitrary 1-dipole contractions, 
that is, combinatorially equivalent graphs (see Definition 4.8) have the same degree. For any 
graph Q we reduce a maximal number of 1-dipoles to obtain a simpler graph with the same 
amplitude as G [86]. 

Combinatorial bubble routing: 

• Designate a root. For a given color i, we pick one of the D-bubbles as a root 7^^-^^ 
bubble. The total number of roots of a graph is 7^^^' = D + 1. 

• Identify i- connectivity graph. We associate to the bubbles i of ^ an i- connectivity graph. 
Its vertices represent the various bubbles ■S^p)- Its lines are the lines of color i in Q. They 

either start and end on the same bubble ^S^^^, known as tadpole lines in the connectivity 

graph, or they do not. A particularly simple way to picture the i-connectivity graph is to 
draw Q with the lines j ^ i much shorter than the lines i. 

• Choose a tree. We choose a rooted tree T* in the z-connectivity graph, such that its root 
is "^(1)- We refer to the rest of the lines of color i as loop lines. 

• Contract. All the iS^ - 1 lines of T are 1-dipoles and we contract them. We end up 
with a connectivity graph with only one vertex corresponding to the root bubble 7^(i)- 
The remaining lines of color i cannot be contracted further (they are tadpole lines in the 
connectivity graph) . The number of the D-bubbles of the other colors is unchanged under 
these contractions. 

• Repeat. We iterate the previous three points for all colors starting with D. The routing 
tree is chosen in the graph obtained after contracting T^^^, . . . ,T^. The number of 
bubbles of colors q > j are constant under contractions of 1-dipoles of color j, hence the 
latter cannot create new 1-dipoles of color q. Reducing a full set of 1-dipoles indexed by 
-D + 1 routing trees T^, . . . , we obtain a graph in which all bubbles are roots. This is 
called a combinatorial core graph. 

Definition 5.1. A combinatorial core graph at order p, denoted Qp^\ is a (I? + l)-colored graph 
with 2p vertices, such that for all colors i, it has a unique D-bubble T^^i)- 
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In fact, note that we already used the combinatorial core graphs three times: first to 
prove (4.2), second to prove Lemma 4.5 and third to prove Proposition 4.3. The important 
feature of combinatorial core graphs is that their degree (and consequently amplitude) admits 
a bound with respect to the number of vertices 2p. By Lemma 4.2: 

i;p 

Combinatorial core equivalence classes: The combinatorial core graph one obtains by the 
above routing procedure is not independent of the routing trees. The same graph leads to 
several combinatorially equivalent core graphs, all at the same order p and all possessing the 
same amplitude. We denote such equivalence by Q^f^ Qp'^^' and call two such core graphs 
combinatorially core equivalent. 

An arbitrary graph Q routes to a unique combinatorial core equivalence class, [^p^^ There 
exists a relation between their respective amplitudes: 

A{g) = (AA)^'"'-^-i^(gW), P = ^ - -D-1). 

Thus, the expansion of the free energy of a (D + l)-colored tensor model can be recast in 
combinatorial core equivalence classes: 

oo 

^(A,A) = E E 

where C^^p^\\,\) counts all the graphs routing (via a combinatorial bubble routing) to the 
equivalence class [^p'^^]- The crucial feature of the expansion in combinatorial core classes is that 
it can be evaluated order by order. One needs to draw all core graphs at order p, compute their 
amplitudes, and divide them into combinatorial equivalence classes. The last step, identifying 
the combinatorial equivalence classes, is potentially difficult, but one needs to deal with this 
problem only a finite number of times to write all terms up to a given order. 

5.5 Topological expansion 

To write a series taking into account not only the amplitude but also the topology of the graphs 
one must rely on topological equivalence rather than combinatorial equivalence. This means that, 
instead of contracting all 1-dipoles, one should contract only 1-dipoles separating spheres [86]. 
Of course it is out of the question to reduce all such 1-dipoles for an arbitrary graph (one would 
need to deal with a very difficult problem for every graph Q). We will modestly only reduce the 
1-dipoles of degree 0. Not only such contractions preserve the topology (from Lemma 4.5), but 
also they lead to a well defined 1/A^ expansion. Unsurprisingly the final series will have some 
degree of redundancy which we will detail later. 
Topological bubble routing: 

• Designate a root. Consider a color i. If there exist D-bubbles of non-zero degree, that is 
cj(;B^^^) > 1, we set all of them as roots Otherwise, we choose one of them as the sole 

root '^(1)- No matter what the scenario, we call 7^^-^^ the principal root and /i > 1 

branch roots. Iterating for all colors, we identify all the roots of G- We denote by T^W the 
number of roots of color i and T^t^^ = 7^[*1 . 



E 
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• Identify i- connectivity graph. As in the combinatorial case. 

• Choose a tree. We choose a rooted tree 7"* in the i-connectivity graph. Consider any branch 
root "^(^^j ^ > 1. It is represented by a vertex in T^*\ We represent by a dashed line, the 

line incident to this vertex and belonging to the unique path in connecting '^(^■j to the 

principal root 'TV^iy All the other lines in 7~* are represented as solid lines. An example is 
given in Fig. 11. 




Figure 11. A tree in the 012 connectivity graph. 



• Contract. All the - 71 ^ sohd hnes in V are 1-dipoles of degree and we contract 
them. We arrive at an i-connectivity graph with vertices corresponding to the roots T^^^)- 
The remaining lines of color i cannot be contracted further. They are either tadpole lines 
or dashed lines. The tadpole lines are not 1-dipoles and the dashed lines now separate two 
roots and so are not 1-dipoles of degree 0. Neither the number nor the topology of the 
bubbles of the other colors is changed under these contractions. 

• Repeat. We iterate for all colors starting with D. The routing tree 7"-' is chosen in the 
graph obtained o/ier contracting T^~^^ , . ■ ■ ■ The degrees (^(7^^^^, q> j are unaffected 
by the contractions of 1-dipoles of color j (by Lemma 4.4), hence the latter cannot create 
new 1-dipoles of degree and color q. Reducing a full set of 1-dipoles of degree indexed 
by D + 1 routing trees , . . . , , we obtain a graph in which all bubbles are roots. This 
is called a topological core graph. 

Definition 5.2. A topological core graph at order p, denoted , is a (D -|- l)-colored graph 
with 2p vertices such that for each color i: 

• either Q has a unique L>-bubble 7^*^-^^ of degree uj{Tl'^^^-^) = 0; 

• or all bubbles 7^(p)> have degree ^(7^*^^^ > 1. 

Just as for combinatorial core graphs, the degree and amplitude of topological core graphs 
admit a bound in the number of vertices. By Lemma 4.2: 

The definition of a topological core graph implies: Ei:p^(^(p)) > - {D + 1), hence: 

u^{gi'^)>[^^^^-iyp + D-ni^^)+p-i>p-i, 

where for the last inequality we used again (4.2). 



36 



R. Gurau and J. P. Ryan 



Topological core equivalence classes: The core graph one obtains by the above routing proce- 
dure is again not independent of the routing trees. The same graph leads to several topologically 
equivalent core graphs, all at the same order p and all possessing the same amplitude. We denote 
such an equivalence by Qp^ Qp^' and call two such graphs topologically core equivalent. 

An arbitrary graph G routes to a unique topological core equivalence class, [Qp 
exists a relation between their respective amplitudes: 



(*)i 



There 



AiG) = (AA) 



P 



2 



Thus, the expansion of the free energy of a (D + l)-colored tensor model can be recast in 
topological core equivalence classes: 



p=l ui>p—l 



E 



[g(*'],a;(g(*')=a; 



(5.2) 



where C^^p ^^(A, A) counts all the graphs routing (via a topological bubble routing) to [Gp^]- The 
classes [Gp^] can again be listed order by order. As before, identifying the topological equivalence 
classes at order p is difficult, but one needs to deal with this problem only a finite number of 
times. 



5.6 Redundancies in the series 

The topological expansion (5.2) groups together only graphs with both the same topology and 
amplitude. However, the scaling with only separates graphs of different amplitudes, as cap- 
tured by the combinatorial expansion. Classifying further the graphs with the same amplitude 
into topological core equivalence classes is a choice. While this leads to a well defined series, which 
one can compute order-by-order, the topological expansion has a certain degree of redundancy: 

• At a given order p, there exist graphs Gp \ Gp^' with the same amplitude but which are 
not topologically core equivalent since they are not topologically equivalent. In formulae: 
^(gW) ^ A{Gi^^') but g^*^ ^ G?' since 6?^*^ oo Although the scahng with N does 
not distinguish two such graphs, we obviously group them in different terms in (5.2). 

• At a given order p, there exist graphs Gp \ Gp^' which are topologically equivalent but 
not topologically core equivalent since they do not have the same amplitude. In formulae: 
G? ~ G?' but Gp^ ?^ G?' since A(af^) ^ A{Gi*'^'). In other words, the same topology 
can appear in several distinct topological core equivalence classes at a given order. 

• For two distinct orders p < q there exist graphs Gp \ Gq^ which are both topologically 
equivalent and have the same amplitude but are not topologically core equivalent since 
they are at different orders. In formulae: Gp^ ~ G^"* and A{G'f') = A{Gf') but Gp'^ ?^ G^"* 
since p ^ q. For instance, if Gq^^ can be obtained from Gp^ by creating 1-dipoles of non-zero 
degree which separate spheres^^. 

• Any topology appears at arbitrary order. For all topological core graphs G? one obtains 
topologically equivalent core graphs (but diminished in power counting) by creating an 
arbitrary number of A;-dipoles of degree for k >2, which separate spheres. 



^^Note that, as the degree of the D-bubbles is additive under 1-dipole creations, q < p + to{Gp) 
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As topology in higher dimensions is a very involved subject, it is to be expected that the 
final expansion is somewhat difficult to handle. Counting all graphs which route to a given 
core graph is relatively straightforward (by performing all possible creations of 1-dipoles of 
degree 0). However, one must pay attention to double counting problems when grouping together 
the graphs routing to a topological core equivalence class. Fortunately, at low orders, the 
topological core equivalence classes have a single representative, and the over counting problem 
is absent. In fact, later in this review, we shall explicitly sum the leading order graphs and 
analyze the critical behavior of the leading order contribution. 

To conclude this topic, let us consider for a brief moment the 3-dimensional case. A bubble 
of a (3 + l)-colored graph Q represents a sphere if and only if = 0. Hence, a 1- 

dipole di separates a sphere if and only if oj{di) = 0. This eliminates most of the redundancies 
of the expansion. 

5.7 First terms of the expansion 

To give us a feel for the expansions describes above, let us calculate the first few terms. Note 
that all combinatorial core graphs are also topological core graphs. The topological core graphs 
up to p = 2 are represented in Fig. 12. In these cases, since they only have one I?-bubble B"^, 
for each color f, they are also combinatorial core graphs. 












Figure 12. Core graphs at p — 1 and p = 2. 

At p = 1, we have a unique topological/combinatorial core graph, denoted g[^^ = Qf\ It 
is the unique topological/combinatorial core graph of degree and thus it represents a sphere. 
For reasons which will become clear below we call this graph the super-melon. As the degree 
of a graph is invariant under arbitrary 1-dipole contractions, the combinatorial class is 
identical with the topological class [Gi^]- 

At p = 2, we have the topological/combinatorial core graphs = with 2<k< [^J . 

This time the combinatorial and topological classes [^2-^^], [^2*1] different. All these core 

graphs are dual to spheres (thus, Q^^ ^2*1;)- 
The amplitudes of the various graphs are: 

The two 1/A^ expansions of the free energy are: 

F(A,A) = C[^i^'](A,A)A^^+ ^ c[^S](A,A)A^^-(^+i)'=+^-'+^ + 0(A^-2), 

k=2 

F(A,A) = C[^i*'](A,A)Ar^+ ^ c[^S](a, A)A^^-(^+i)'=+'='+^ + 0(Ar^-(i^2^. ^^^3^ 

k=2 
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Note that the combinatorial expansion has a better a priori bound on the remainder terms than 
the topological expansion. This should come as no surprise, as the topological classes [^2*1;] 

contain fewer graphs than the combinatorial classes [^2^^^,]- At first sight, the bound on the 
remainder terms is not very tight. The first explicit correction with respect to the dominant 
behavior in (5.3) comes from ^2,2- To ensure that all terms with amplitude larger or equal 
than ^2,2 have been taken into account, one needs to check combinatorial core graphs up to 
Pmax,(c) = -D — 1 in the combinatorial expansion and topological core graphs up to Pmax,(t) = 
1 + ~ ^) ™ topological expansion. 

6 Embedded matrix models 

We return here to the graph factorizations we developed in Section 4. Our earlier disposition 
on graph factorization was intended to highlight the relationship between factors of labeled 
complete graphs JCd+i and subsets of {D + l)-colored graphs G. We saw that 1- and 2-factors 
translated into patches and jackets, respectively. In particular, a factorization of /C^+i in terms 
of 1- and 2-factors determined a P J- factorization of Q. 

The colored version of ICd+i makes a reappearance when we examine the tensor model (5.1). 
Irrespective of the index structure, each interaction term clearly generates a vertex with D+1 
distinctly colored incident lines. The boundary of such a graph component is a colored /Cd+i- 
Thus, at the level of the action, we have the same picture as earlier: a single {D + l)-valent 
vertex and its boundary graph. 

The colored tensors are associated to the edges of the interior and thus the vertices of the 
boundary graph. Our choice of index structure rests on which substructures we would like to 
weight in our tensor model. In the i.i.d. action, each tensor has D indices and a given tensor is 
contracted with each of the other D tensors using one of these indices. Given that the tensors 
are associated to the vertices of the boundary graph, these contractions weight the edges in the 
boundary graph. As a check, from the counting above, we have D(D + l)/2 pairs of contracted 
indices, which coincides with the number of edges in ICd+i- 

Matrix models and jackets: A 2-factor of /C_d+i is a (-D + l)-cycle in the vertex set of /Cd+i 
and thus picks out + 1 edges. Translating this into tensor language, this cycle picks out D + 1 
pairs of contracted indices in the interaction term and each tensor contains one element of two 
pairs. Neglecting all other indices for the moment, this has the character of a (D + l)-colored 
matrix model. To illustrate this, let us examine the 2-factor (0, 2, 1, 3, 4, . . . , D). It picks out 
the edges (02), (21), (13), (34), and so on up to (DO). Then, the interaction term in the action 
takes the form: 

A ^ 

n i=A 

^ Tr^(<AVW0^---<A^)+c.c., 



^Z)(D-l)/4 



where we left all other indices implicit and introduced some new notation Trj- meaning the 
trace over the indices associated to the edges picked out by the 2-factor. As we have written 
it, this subset of the tensor model is a matrix model in its own right. As a result, it generates 
a Riemann surface none other than the jacket J determined by the 2-factor. This is a general 
rule: jackets are weighted by matrix models embedded inside the tensor model. 

Vector models and patches: For odd D, a 1-factor is a pairwise association of vertices of JCd+i 
and thus picks out {D + l)/2 edges. Translating this into tensor language, the 1-factor couples 
the D + 1 tensors into {D + l)/2 pairs. Neglecting all other indices for the moment, this has 
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the character of a set of coupled vector models. Consider the 1-factor specifying the edges 
{{a + {D + l)/2) : i G Z(£)_(_x)/2}- Then, the interaction term takes the form: 

D-l D+l 
n ^ 2 2 

(/ -0^) ••• (^V .^D) +C.C. 

Trp +C.C., 

where in the second line we see its manifestation as coupled vectors and in the third line, we 
have introduced another piece of notation Tr-p meaning the trace over the indices associated 
to the edges picked out by the 1-factor. These vectors are, once again, embedded inside the 
tensor structure. This vector structure does not generate a Riemann surface this time, but it 
still generates the faces of the patch V determined by the 1-factor in question. This is another 
general rule: patches are weighted by vector models embedded inside the tensor model. 

We now possess enough tools to translate the concept of PJ-factorization into the tensor 
model setting. As we have stated, a PJ-factorization is determined by a factorization of edge 
set of JCd+1 into 2-factors (and, for odd D, a single 1-factor). In our tensor model, contracted 
pairs of indices correspond to edges in /C/5+1, so that we can factorize them into 1- and 2-factors. 
Thus, for even D, we find that indices of the interaction term are exhausted by D/2 2-factors, 
while for odd D, they are exhausted by {D — l)/2 2-factors and a 1-factor. 

This factorization allows us to rewrite the action as: 

S{cp, 0) = E Trim + ^^(ii)/4 Trp^ (U A + ^z>(D-i)/4 ^^^J (U A ' (6-1) 

i=0 \i=0 / \?=0 / 

where Tr-pj- denotes the matrix trace over the indices associated to the disjoint 2-factors and 
the scalar product over the indices associated to the possible 1-factor. Once a particular PJ- 
factorization is chosen, this partition of degrees of freedom among the jackets is preserved by 
amplitudes. Remember that: 

A{g) = (AAfiV^^P^"'^^^ = (AAf A^^(^^e)^ 

using equation (4.6). This amplitude is a remarkably simple generalization of the 2-dimensional 
case. We note that each jacket in the factorization contributes a factor N'^^^\ exactly what one 
would expect from the matrix model associated to Riemann surfaces. For odd dimensions, there 
is also a patch which contributes a factor xi'P)- While such coupled vector models are not so 
common in the literature, the amplitude xi'P) is appropriate to such a model. 

Remark 6.1. In Section 5, we have developed a 1/A^-expansion such that the sphere dominates 
in the large N limit (due to the fact that it is the only topology that contributes at w(^) = 0). 
In other words, all the jackets have zero genus. But this reasoning was formulated treating all 
jackets in a rather egalitarian manner. Now that we have chosen a particular factorization, we 
see that the dominant configurations correspond to having only spherical jackets and in certain 
cases, a patch such that xi'P) is maximal. In using a particular factorization, we have chosen 
to highlight certain surfaces. One could imagine making this distinction more pronounced, for 
example, by weighting the different surfaces by different scaling parameters. In such a sce- 
nario, the tensors 0* are no longer maps with the domain (Z^r)^, but rather with the domain 
XjZatj, X ZjVp, where J'i and V are respectively the jackets and the (possible) patch in a PJ- 
factorization. We could then scale each of the jackets independently. In this case, no longer do 
the configurations with w(Q) = necessarily dominate. We could for example make just Nj 



40 



R. Gurau and J. P. Ryan 



for just one jacket large, keeping the other relatively small. In this case, all configurations with 
one spherical jacket would be equally dominant. In Proposition 4.3, we showed that all these 
configurations were still spheres, at least in D = 3. Thus, these generalizations are of interest as 
they are likely to exhibit different critical behavior to the simplest i.i.d. model (see Section 7). 

Remark 6.2. The identification of embedded matrix and vector models [119, 120] provides 
a first glimpse of an currently active research avenue in that it hints at how one might adapt 
the plethora of powerful multi-matrix model techniques to the higher-dimensional setting [53]. 
Although it is unlikely that these are utilizable in the most generic cases, there is tantalizing 
evidence that they can be used to solve models exactly in specific regimes [120]. 

7 Critical behavior 

We are finally ready to examine the critical behavior of tensor models. We shall see in this 
section that the colored tensor models undergo a phase transition to a continuum theory of 
large spaces. The topology of the space thus obtained is spherical^'^ The leading order in 
the 1/A^ expansion consists of a subclass of colored triangulations of the i^-sphere. We will now 
address the combinatorics of the dominant triangulations. As usual, we include an abridged 
version of the forthcoming developments here. 

The dominant order: Our first task is to identify the combinatorics of the leading order graphs 
(which are all spherical). It emerges that all graphs have a characteristic 'melonic' structure. 
The leading order graph with the lowest number of vertices is the 'super-melon' presented in 
Fig. 14(a). All other leading order graphs are obtained by repeatedly inserting D-bubbles with 
two vertices (i.e. two vertices connected by D lines) arbitrarily on lines. As we explain in a short 
prologue to this subsection, it is possible and more convenient to deal explicitly with 2-point 
graphs and translate to vacuum graphs at the end. 

From melons to trees: The melonic graphs are in one-to-one correspondence with the so- 
called colored rooted {D + l)-ary trees. We present their definition and a detailed map between 
the two combinatorial species. 

Resumming the dominant series: Having recast the dominant order as a series indexed by 
a species of trees (and possessing a finite radius of convergence), we present no less than three 
ways to perform the arithmetic: a direct method; a method using the {D + l)-ary tree measure 
from Section 2; and a method using Cayley's theorem. 

Critical exponent: We extract the critical value of the coupling constant and critical exponent 
associated to the free energy (the 'string' susceptibility exponent 7). 

Melons as branched polymers: We outline the interpretation of the (D-l-l)-colored i.i.d. tensor 
model as the generator of dynamical triangulations. This allows us to visualize the problem 
from a more geometrical and hopefully physical standpoint. Translating our dominant melonic 
graphs into this formalism, we find that they display a remarkable resemblance to the sector of 
dynamical triangulations known as branched polymers. They both have large positive curvature 
and moreover their entropy exponents match. However, the spectral and Hausdorff dimensions 
of melonic graphs are still to be determined before concluding on the precise relation with 
branched polymers. 

Large volume limit: We conclude this section with a description of the large volume limit. 
This is obtained by sending the length (a parameter introduced in the map to the dynamical 
triangulation picture) of the individual simplices to zero and the coupling constant of our tensor 
model to its critical value. We highlight the subtlety inherent between taking this limit and the 
continuum limit. For the latter we must also keep the renormalized cosmological constant fixed. 
We finish with some interesting current research directions. 



Its precise geometry, however, is still to be understood. 
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To get to the melon rolling, so to speak, the free energy and the 2-point function of the 
i.i.d. model are related by a Schwinger-Dyson (SD) equation. (This is in fact true for a generic 
tensor model.) It turns out that in order to count the leading order graphs, it is easier to count 
the 2-point graphs and recover the free energy using this equation. Due to the conservation of 
the indices along the strands, the 2-point function is necessarily connected and has the index 
structure: 



i^k^k)c = ^ri.,n.G{X,X), (7.1) 

where 6^^^^^^ denotes Uk^i ^n^^^n,, ■ 

Proposition 7.1. The full connected 2-point function is: 

G(A,A) = 1 + XdxE{X,X), 

where E = N^^F is the free energy per degree of freedom of the model. 

The proof of this, by standard techniques can be found for instance in [33]. In words, this 
SD equation tells us that vacuum graphs can be obtained from 2-point graphs by reconnecting 
the two external edges (and conversely, cutting any edge in a vacuum graph generates a 2-point 
graph). 



7.1 The dominant order: the world of melons 

We present below the precise construction of the leading order melonic graphs. As the construc- 
tion is somewhat convoluted it is useful to first state the end result. Call S the self energy of the 
model (that is the 1 particle irreducible amputated two point function). At leading (melonic) 
order Smcionic factors into the convolution of D connected two point functions Gmcionici one for 
each color. This is represented in Fig. 13. We will now establish this result. 



G 



rnclonic 




melonic ) 



^melonic 



Figure 13. The IPI 2-point function and the full connected 2-point function for melonic graphs. 

Our construction relies on eliminating D-bubbles with two vertices. A D-bubble with two 
vertices is a D-dipole^^. Thus, by Lemma 4.4, the degree and topology of a graph Q are 

identical to those of the graph Q/B^^^^ obtained by replacing S^^-j with a line of color i (see 
Fig. 14(b)) 

For D > 3, it can be shown [33] that a leading order 2-point graph must possess a D-bubble 
with exactly two vertices. First one proves: 

Lemma 7.1. If D > 3 and Q is a vacuum graph with degree 0, then Q has a face with exactly 
two vertices. 



This is the case unless it is a bubble of the super-melon graph depicted in the Fig. 14(a). 
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i + l 



D i-1 
(a) The super-melon graph. (b) A D-bubble with two vertices. 

Figure 14. The super-melon and D-bubble elimination. 

Proof. Since the graph Q is of degree 0, from equation (4.3) we can ascertain that its face set 
has cardinality \Tg\ = ^ ^^^ p + . Let us denote by the number of faces with 2s 

vertices (every face must have an even number of vertices). Then: 

iJ'gM + \^g,2\ + Yl = + D. (7.2) 

s>3 

Let 2|?|'^^ be the number of vertices of the p-th face with colors ij. On the one hand, we count 
the total number of vertices by summing the numbers of vertices per face. On the other hand, 
each vertex contributes to D{D + l)/2 faces. This leads us to: 



— ^ — p- 



We solve (7.2) for the number |-Fg.2| of faces with four vertices, and insert the result in (7.3) to 
get: 



\Fg,,\ =2D + Y,is - 2)\Tg,s\ + ^^^^i 



-P- 

z 

s>3 

Notice that on the right hand side, the first two terms yield a strictly positive contribution for 
any D >2, whereas the third term changes sign at D = 3. Thus, we conclude that: 

|J"g,i|>l if L>>3. ■ 

Using Lemma 7.1 and the fact that all jackets of a leading order graph are planar, one can 
prove [33]: 

Proposition 7.2. If D >3 and G is a leading order 2-point graph, then it contains a D-bubble 
with exactly two vertices. 

This result reveals some properties of leading order vacuum graphs. Consider a leading order 
vacuum graph Q. It is either the super- melon graph of Fig. 14(b) or it contains a D-bubble with 
exactly two vertices which can be eliminated. This elimination procedure preserves colorability, 
degree and topology and so we obtain a leading order vacuum graph with two less vertices. Thus, 
the new graph satisfies the same conditions: either it is the super-melon graph or it possesses 
a D-bubble with two vertices. In the latter case, we iterate the procedure. It follows that the 
leading order vacuum graphs reduce to the super-melon graph after sequentially eliminating 
their D-bubbles. When counting the leading order graphs it is more useful, however, to take the 
reversed point of view. One starts with the super-melon graph and inserts D-bubbles with two 
vertices on its lines. 

Now let us move back to dealing exclusively with leading order connected amputated 2-point 
graphs. 
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Definition 7.1. A melon graph 7W is a leading order 1-particle irreducible (IPI) amputated 
2-point graph. 

Not all leading order connected amputated 2-point graphs are IPI (see for instance Fig. 17), 
but one may build a generic leading order 2-point graph from those possessing this property. 
Thus, it is possible to introduce the following class of graphs. 

Definition 7.2. A melonic graph ^ is a graph in which all IPI amputated 2-point subgraphs 
are melon graphs. 

All leading order amputated 2-point graphs are melonic graphs (although not necessarily 
melon graphs, as they are not necessarily IPI), as is the case for the example drawn in Fig. 17. 
Additionally, we shall denote the set of melons of G hy Sg. 

The intuitive picture is that a melonic graph itself comprises of melons within melons. The 
Z?-bubble with only two vertices is obviously the smallest melon, while for the leading order IPI 
amputated 2-point graph, the largest melon is the graph itself. 

At leading order, the connected 2-point function and its related free energy take the form: 

oo oo 
Cmelonic = G'p(AA)^, -E-melonic = ^^-E^p(AA)^, 

p=0 p=0 

where Gp (resp. Ep) is the number of leading order 2-point (resp. vacuum) (D + l)-colored 
melonic graphs with 2p vertices. We want to find a closed form for Gmcionic (and hence -Emcionic) 
by resumming the above series. 

Direct solution: We now extract some juice from the world of melons. Let Smeionic be the 
melonic IPI 2-point function (with external edges of color D). Any contribution to Smdonic can 
be obtained by arbitrary insertions of melonic graphs on the interior edges of a D-bubble. (We 
do not allow insertions on the active exterior edge of the root melon, since this would cause 
the graph to lose its IPI property). Thus, Smeionic necessarily has the structure represented 
in Fig. 13. Taking into account (5.1) and (7.1), it can be written as a function of the melonic 
2-point function Gmelonic: 

~ ^ ^ X/ n ( ^melonic )'^ = (AA) 5^^^^ (Gmelonic)^ • 



JVD(Z)-l)/2 



One gets a closed equation for Gmcionic by recalling that the sum of the geometric series of IPI 
amputated 2-point functions yields the connected 2-point function: Gmeionic = (1 ~ 5]jnciomc)~^- 
Hence, 

Cmelonic = 1 + (AA) (Gmclonic)'^^^ • 

The condition that at leading order in the large limit the IPI two point functions factors 
into contributions from the connected two point functions is at the core of all the studies of the 
critical behavior of tensor models performed so far. 

This equation is well-known in the literature ([74, Exercise 2.7.1], [75, p. 200], [126, Propo- 
sition 6.2.2]) appearing in various problems of enumeration [92]. Similar equations appear in 
the study of branched polymers. Consequently, the leading order in 1/A^ shares the same 
combinatorics with the branched polymers, hence many of their statistical properties. Such 
configurations have been observed in numerical studies of dynamical triangulations. However, 
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as already mentioned, there is a major point of divergence between branched polymers and our 
melonic graphs. Namely, the intrinsic Hausdorff and spectral dimensions of the melons are yet 
to be computed, and as we will explain below, one has no a priori reason to expect that they 
will equal those of the branched polymers. 

The solution of equation (7.4) which goes to 1 when (AA) goes to zero can be written as 
a power series in (AA) with coefficients the {D + 1)-Catalan numbers. 

Proposition 7.3. The melonic 2-point function admits the following expansion: 



G„,eionic(A,A) = J^C7(^+i)(AA)^ 

p=0 



wtth c(^+i) = — ^- 

P {D + l)p+l 



{D + l)p+l 
p 



the (D + \)-Catalan numbers. 

The usual Catalan numbers correspond to the case D = 1. Note that the case D = 3 can be 
solved explicitly: 



G 



molonic|_D=3 



(l + 4^;)i/4_ [2-(l + 4t;)V2]'^' 



where 



9^ 
2V3 



1 + A 1 



28 
33^ 



1/3 



+ 1 



28 
33^ 



1/3- 



(7.5) 



and 



ff = AA. 



We shall see below that the leading order connected 2-point graphs are in one-to-one corre- 
spondence with colored rooted {D + l)-ary trees [109]. 



7.2 From melons to trees 

We calculate directly the coefficients Gp. To do so, we pick a fixed color, D say, for the pair of 
external legs and count the number of melonic graphs at a given order. Remember that once we 
list the dominant graphs at a given order, we can generate those at the next order by arbitrarily 
inserting D-bubbles with two vertices. In fact, the lowest orders in p can be evaluated by a direct 
counting of Wick contractions: Go = Gi = 1, G2 = {D + 1) and so on, but since we need all 
orders, we provide a systemic treatment. 
Order \\: 

• The lowest order graph consists in exactly one Z?-bubble with two vertices (and external 
lines of color D). There is only one Wick contraction leading to this graph, hence Gi = 1. 




Figure 15. The first order melonic graph and its corresponding rooted tree. 

• We represent this graph by the rooted tree with one central vertex decorated with D + \ 
leaves and one root. A leaf is a vertex with only one incident edge. These D + \ leaves 
correspond to all the edges touching the positive vertex. The leaves inherit the colors. 
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The root corresponds to the external edge of color D touching the negative vertex. The 
vertex to which the root belong is called the root vertex. We consider all edges incident 
at the positive vertex to be active, while the external edge incident at the negative vertex 
is inactive. The leafs of the tree inherit this (in)activity. See Fig. 15 for an illustration, 
where the inactive line and leaf are represented as dashed lines. 

Order {XXf: 

• At second order, D + 1 graphs contribute. They arise from inserting a -D-bubble with two 
vertices on any of the D + 1 active lines of the first order graph. Say, we insert the new 
bubble on the active edge of color j. With respect to the new Z)-bubble, all edges incident 
at its positive vertex are deemed active, while the exterior edge (of color j) incident at 
its negative vertex is deemed inactive. Each graph has a combinatorial weight but 
it is produced by 2!^ different Wick contractions. These correspond to the independent 
relabeling of the positive and negative vertices. Thus, the overall factor is 1. 




D-1 D 



Figure 16. A second order graph and its corresponding tree. 

• Once again, let us examine the case where we insert the new bubble on the active edge of 
color j. This graph corresponds to a tree obtained from the first order tree by connecting 
its leaf of color j to a new {D + 2)-valent vertex. This new vertex has D + 1 leaves, one of 
each color. The remaining line, known as a tree line, is the line of color j joining the two 
(D + 2)-valent vertices. The leaves correspond to the active lines of the new bubble and 
hence are themselves active. We count D + 1 distinct trees, hence G2 = D + 1. See Fig. 16 
for the case j = 0. We shall introduce a canonical labeling of colored rooted {D + l)-ary 
trees in Section 8.1. With regards to this labeling, each of these trees is a {( ), (j)} tree. 

Order (AA)p+^- 

• We obtain the graphs at order p + 1 by inserting a Z?-bubble with two vertices on any of 
the active lines of a graph at order p. Once again, with respect to the new bubble, all 
edges incident to its positive vertex are deemed active, while the exterior edge incident to 
its negative vertex is deemed inactive. 

• In terms of the trees, we represent this insertion by connecting a (Z? + 2)-valent vertex, 
with D + 1 active leaves, to one of the active leaves of a tree at order p. The new tree 
line inherits the color of the active edge on which we inserted the L>-bubble. Each of these 
graphs has a combinatorial weight (^p_^iy2 ; while is produced by (p+1)!^ Wick contractions 
(corresponding to the relabeling of of positive and negative vertices). Thus, the overall 
combinatorial weight is 1. 

To recapitulate, at a given order (AA)^, we obtain contributions (with combinatorial weight 1) 
from all rooted colored {D + l)-ary trees with p vertices (and hence with Dp + 1 leaves). We 
remark that the tree associated to a graph is a colored version of a Gallavotti-Nicolo (GN) 
tree [72]. 
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A tree is associated to a partial ordering. In our case this is an ordering, denoted >, of the 
set Sg of melons of the graph. For two melons Mi,Ai2 G Sg: 



Ml > M2 if 



'either Mi D M2, 

or 3M^p)£Sg:gDMiO ( (jM^p) ] D M2, 



where we have defined a U operation on Sg such that M1OM2 denotes the connection oi Mi 
and M2 using the active external edge of Mi and the inactive external edge of M2- In addition, 
> is transitive. 

We may depict this partial ordering very neatly using the tree graph developed earlier. First, 
not that we can identify each melon by the positive vertex at which its active external edge is 
incident. Since the positive vertices are in one-to-one correspondence with D + 2-valent vertices 
of the tree, we may label each such vertex of the tree graph by the associated melon. Then, the 
tree graph captures the partial ordering explicitly. The root melon is the melon mapped to the 
root vertex of the tree. For a melon M2, if another melon Mi lies on its unique path to the 
root melon, then Mi > M2- Mi is an ancestor of M2, while conversely M2 is a descendant 
of Ml. 

An example in D = 3 is given in Fig. 17. The root melon is Mi. Note that M3 D 
M4,M5,Mq,M7, hence it is their ancestor. Also, since MsOMsOMiq C G and M9 C Mio, 
the melon M3 is the ancestor of Mi, M5, Me, M7, Ms, M9, Miq. 




Figure 17. A melonic graph and its associated colored GN rooted tree. 

The leading order connected vacuum graphs (also called melonic) are obtained by reconnect- 
ing the two exterior edges of a melonic 2-point graph. Their amplitude is A^^. If a graph is 
a (-D + l)-colored melonic graph, all its D-bubbles are melonic graphs with D colors. This is easy 
to see, as the reduction oi a D bubble with two vertices represents the reduction of a (L* — 1)- 
bubble with two vertices for all D-bubbles which contain these two vertices. When reducing the 
graph to its root melon, one by one all its -D-bubbles reduce to D-bubbles with two vertices. 

Moreover, the D-ary trees of the -D-bubbles are trivially obtained from the (Z) + l)-ary tree of 
the graph by deleting all lines (and leaves) of color D. We shall need below the following obvious 
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fact: given a melonic graph and one of its D-bubbles (say all the lines of color D connecting 

on it either separate it from a different D-bubble (and are tree lines in the associated colored 
rooted tree) or they connect the two external points of a IPI amputated 2-point subgraph with 
D — 1 colors of (i.e. they connect the two external points of a melon in ^ji^): in which case 
they are leaves of the associated tree. 

We shall display below two additional techniques to resum the melonic graphs using this 
correspondence with trees, since they are important tools in their own right. 

Counting using the D + 1-ary tree measure: From Section 7.2, we conclude that Gmeiomc('^-^) 
is the generating function of rooted {D + l)-ary trees. The equation (7.4) encodes their pro- 
liferation. Its interpretation in term of trees is straightforward, and depicted in the Fig. 18. 
The leading order connected 2-point function expands in colored rooted trees with p (unlabeled) 
D + 2-valent vertices, Dp + 1 active leaves, and one inactive leaf taken as the root. The result 
is recovered using Example 2.7. 



O 



1 



+ AA 




Figure 18. The melonic equation represented from the point of view of the colored rooted trees. 

Counting using Cayley's theorem: Cayley's theorem states that the number ordinary trees 
on n labeled vertices, such that the vertex i is dj-valent, is: (n — 2)!/ J3j((ij — 1)!. We shall work 
towards applying this result bit by bit. First let us denote by Tp, the number of ordinary trees 
with p labeled {D + 2)-valent vertices and Dp + 2 labeled 1-valent vertices. Using the above 
formula, one has: 

_ [{D + l)p]\ 

There are two differences between ordinary trees and our (L>+l)-ary trees: i) the vertex-labeling 
and ii) the line-coloring. We first start by coloring the tree. The color of the inactive leaf of the 
root vertex is fixed by the external color of the 2-point function. We have {D + colorings of 
the lines of the root vertex. Any {D + 2)-valent vertex will support {D + \)\ colorings of its lines 
(as the color of the line connecting it towards the root is fixed by its ancestor) . Thus we pick 
up a factor [[D + 1)!]^ for the colorings of the lines. We then remove the vertex labels. This 
brings a factor ^ for the {D + 2)-valent vertices and a factor ^j^^j^^y^ for the active leaves. We 
conclude that: 

_ [{D + iy]' _ 1 1 _ (^+1) 

" p\{Dp + 1)! " p\ {Dp + 1)! + ~ ■ 

Note that there is one subtle point in using Cayley's theorem: one can not first relabel the 
vertices and then color the lines. Indeed, in order to identify the allowed colorings of the lines 
one must have distinguished vertices (otherwise one can not properly count for instance the 
colorings of leaves touching the same vertex) . 

7.3 Critical exponent 

Of course the whole point of this counting was to use the outcome to categorize some critical 
behavior of the class of models in question. Applying Stirling's formula to the expansion in 
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Proposition 7.3, we get the large order behavior. 

Proposition 7.4. The number of melonic graphs with 2p vertices behaves for large p like: 
with 

(Z) + 1)^+1' """^ 1)3/2 ■ 

In particular, Gp is exponentially bounded by e"^^^^", with Inge < 0. 

The D-dependent constant gc is the critical value of the coupling g = XX. Indeed, it is well- 
known that a series with coefficients of the form gc^P~" behaves in the neighborhood of gc like 
((7c — g)°'~^. Hence, the most singular part of the melonic 2-point function is: 



^ 1/2 

9c- g' 



Cmelonic,sing ~ ^ 

\ 9c 

for some constant K}^ The free energy of melonic graphs is obtained using Proposition 7.1: 

(g \ ^ O'mcloilic ^ 

1 ) Tmelonic — 7T' 

9c J ^ 

The exponent 7meionic is known as the susceptibility, or entropy exponent^^. 

Alternatively, one can derive the critical exponent from equation (7.4), noting that its solution 
may be written as: 

a = G~^ - C~^~^ 

" melonic melonic ' 

and that in the neighborhood of the critical point: gjj-^^ — "(^melonic) ~ 0' ^^sX the coupling 

expands as: g - gc = (Gmcionic - Gj^eionic)^- 

This critical behavior is a key ingredient in the provision of a continuum limit for random 
colored melonic triangulations. But, a priori, the geometric interpretation of the continuum 
limit depends on the details of the model under consideration. In the following, we use the 
natural interpretation of the i.i.d. model as a generator of dynamical triangulations. 

7.4 Melons and branched polymers 

One can interpret the i.i.d. model quite naturally as a model of dynamical triangulations (DT) 
[2, 3, 5, 11, 48]^^. In the remainder of this section, we switch to notations more familiar in 
the DT literature. A melonic graph is dual to a colored triangulation of the D-sphere. We 
denote the number of fc-simplices (i.e. the number of (D — A;)-bubbles) by N^- In this case, the 
amplitude of a graph G becomes: 

A{g) = e>'0-2ND-2-KoND ^ (7,g) 

where K£)_2 and kd are: 

KD-2 = lniV, and kd = ]^{^D{D - l)\nN -\n{g) 



^^This can be checked explicitly from (7.5) in the D = 3 case. 

^^Once an exponential bound is found for the proliferation of a species, the entropy exponent characterizes the 
polynomial part of the number of such objects. 

^ '^Tensor models were indeed originally proposed in this context [7, 124]. 
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Once again, g = AA. The argument of the exponential takes the form of the Regge action (dis- 
crete form of the Einstein-Hilbert action for general relativity) on a triangulation with regular 
L>-simplices of length, say, a. Indeed, on a Regge discretization, the curvature is concentrated 
around the {D — 2)-simplices and measured by the deficit angle 5{aD-2) (2vr minus the sum of 
the dihedral angles hinged on the {D — 2)-simplex 00-2) and the total volume is measured by 
the cosmological term. The Regge action is: 

^Regge = A^VoI(o-d) - ^ Vo1(c7D-2)(^(o-D-2) • (7.7) 

CO crD-2 



The volume of a regular fc-simplex is: vol(o"fc) = xTy^^' ^^"-^ (7-7) takes the form (7.6) on 
a regular triangulation if we identify N and g in terms of the (bare) dimensionful parame- 
ters G, A and the length a via: 

\ng = vol(cjD-2) (^{D -1)-{D + 1) arccos ^ - 2A voI(ctd) = -2a^ A, 

lniV=^^°^. (7.8) 

Notice that in two dimensions, Ing = — 2A vo1(iT£))- 

The large limit corresponds to G — )• 0. Since g is kept finite, A becomes large, positive, 
and scales like A ~ l/(a^G). The melonic triangulations have degree 0, hence their amplitudes 
take the form: 



where we have explicitly split the contribution of the Einstein-Hilbert term and of the cosmo- 
logical term. The scalar curvature is the argument of the first exponential (up to 1/G and 
irrelevant constants). It grows linearly with the number of D-simplices and is positive (since 
tt{D — 1) — (D + 1) arccos ;^ > for D > 3). The entropy, that is the logarithm of the number 
of triangulations with a fixed volume (fixed number of D-simplices) , grows linearly with p for 
melonic graphs, by Proposition 7.4. 

For the melonic family, the numbers of D-simplices and (D — 2)-simplices are not independent. 
They are related by: 

Nd-2 = ^^ Nd + D. 

Hence, the large limit of the colored tensor model projects dynamical triangulations to curves 
(parametrized hy g = AA): 

D{D-l) 1 1 
^ = -- In^ > -- In^c, 

where gc is given in the Proposition 7.4. This gives the critical curve: 
Cl , D{D-1) 1 

Kd{kD-2) = ^ KD^2 - - Inffc- 

The family of melonic triangulations shares a number of similarities with the branched poly- 
mers (BP) phase of dynamical triangulations [3, 4, 6, 11, 48]. Indeed, BP are known to dominate 
the regime of large positive curvature, > i^i){hd-2) for sufficiently large values of kd-2- In 
Fig. 19 we present the phase portrait of dynamical triangulations. The continuum phase is 
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Figure 19. The phase portrait of dynamical triangulations. 

reached for k^(k£)_2)- The BP phase corresponds to > f^%^2- dynamical tri- 

angulations undergo a (first order) phase transition to the crumpled polymer phase. Since 
kd-2 = In A^, one would expect that the large limit of tensor models yields the BP phase of 
dynamical triangulations. 

The case D = 2 is special. Although melonic graphs appear in the large A^ limit, the critical 
behavior is driven by the whole planar sector [28] (and not only by the melonic graphs). Thus 
Tstring = — |- Only beyond the c = 1 barrier the dominant family is reduced to BP, with 
a susceptibility exponent 7bp = 2- higher dimensions, the large A^ limit yields directly 
a phase of large positive curvature with critical exponent 7meions = 7bp = ^, hence with the 
statistical properties of BP. 

The last comparison between melonic triangulations and BP is geometric. Indeed, it has been 
observed in DT that in the BP phase the number of vertices and the number of D-simplices 
grow proportionally A'^q — No/D (see for example [48]). The natural geometric interpretation 
is the following. Start with a triangulation of the D-sphere and iterate random sequences of 
(!—)•(£' + 1)) moves. Such moves add a vertex to a L>-simplex and split it into {D + 1) new D- 
simplices. This way, the triangulation gets a new vertex together with D additional D-simplices. 
Hence, asymptotically, A'o — Nd/D. Notice that in the crumpled polymer phase of DT, one has 
in contrast A^q ~ ^ -^D- Melonic graphs respect: 

hence behave like BP. Irrespective of the dimension, the coefficient 1/2 arises in the first term 
on the right hand side of the above equation. Its occurrence has the following explanation. 
Melons are not obtained just through (1 — s- {D + 1)) moves^*^. Instead, one adds a D-bubble 
with only two vertices to some initial melon by first performing a (1 — t- (Z) + 1)) move, which is 
immediately followed by a (D — )• 2) move^^. Such a sequence adds one vertex to the triangulation 
together with only two D-simplices. Clearly melons provide a better balance between Nq and No 
(maximize Nq for a fixed N^)) than the usual set of BP. It may be surprising that they have 
never been considered as such in DT. We think that this is because the creation of a D-bubble 
with two vertices needs two successive moves at the same place, which is unlikely to happen in 
Monte-Carlo algorithms. 

The main reason we refrain for now from stating that the melonic graphs are BP it the 
following. The BP phase is characterized by a Hausdorff dimension dn = 2. Obviously, one 
could get this value for melonic graphs by using the natural distance on the colored trees. 
However, the vertices of the trees correspond to melons, which are non-local from the point of 
view of the triangulation. One should instead use a different notion of distance (graph distance 
either for the triangulation or for the melonic graph), and not the tree distance in the abstract 
colored tree. 



'This move alone does not respect the colors. 



^^This second move restores the colors spoiled by the first one. 
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7.5 Large volume limit 

The continuum, large volume limit is obtained, like in matrix models, by simultaneously sending a 
to zero and g = XX to its critical value. When g <^ gc, triangulations with a large number of 
D-simplices get exponentially suppressed by {-^)^°^'^ and triangulations with a small number of 
D-simplices dominate the free energy i?meions- When g ^ gc, the summability of -^melons is lost 
due to the triangulations with an infinite (i.e. large) number of simplices, whose individual sizes 
scales to zero. Notice however that in this large volume, continuum limit the coupling hd-2 is 
sent to infinity (and thus G to zero) instead of tuning it to some critical value. This is exactly 
like the standard large N limit of matrix models for pure two-dimensional gravity and here, 
it is simply due to the fact that kd-2 = InA^. To access analytically a sector of finite K£)_2 
would require some kind of double scaling limit [43, 59, 77] which remains still to be uncovered 
in tensor models. 

The average volume is: 

(Vol) = (Nd) = 2a'^A0A InE'meions = "^dc 1- less singular terms, 

gc- g 

and must be kept finite in the continuum limit^^. In the large N limit, we are left with a single 
dimensionful parameter, A in (7.8), which should be renormalized when reaching the critical 
point. The average volume is written in terms of A as: 

2Ar' ^ -Inc/ V 5c / 

In order to obtain a continuum phase with finite physical volume, we need to take the continuum 
limit: 

a^O, g ^ gc, with Ar ~ a~-^(5rc - 9) fixed. 

Notice that this renormalization of A has already appeared in the DT literature, for example [3], 
where it was given the equivalent form: 

o'^Ar = kd - k'o{kd-2)- 

An important issue concerning the critical behavior we have observed is to understand its 
universality, i.e. the relevance of the details of the microscopic model. In matrix models for 
2D gravity, one can describe the surfaces using either triangulations or quadrangulations. Ac- 
tually, any interaction with a positive coupling falls in that same (pure gravity) universality 
class [53, 103]. We will perform the first steps in the study of universality of colored models 
in arbitrary dimensions in the sequel. Establishing the full (multi-) critical behavior of colored 
tensor models is one of the active research directions in the field. 

An interesting feature revealed in our analysis so far is a universality with respect to the 
dimension D. Although the value of the critical coupling depends on D (and goes to zero 
when D goes to infinity), the susceptibility exponent does not. This is in agreement with 
the almost systematic appearance of a branched polymer phase in dynamical triangulations. 
Although we have restricted our analysis to the i.i.d. model, the same melonic graphs dominate 
more involved models like the colored Boulatov-Ooguri model [90]. In fact, as melonic graphs 
have the maximal number of faces at fixed number of vertices, they are likely to generically 
dominate in tensor models, irrespective of the details of the covariance (propagator). Since 



(Vol)~-^(^) 
2A V 5c / 



Like in the case of matrix models, one needs in fact to look at geometries with at least a boundary component 
in order to truly get an infinite volume. 
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we have mapped melons to trees, it is appealing to think that they inherit all the statistical 
properties of trees. However, the mapping does not preserve the locality of the triangulation, 
so that more investigations are needed in order to conclude whether the physical correlations 
defined on the triangulation agree with those of branched polymers. In particular, it would be 
interesting to test Fisher's scaling relation directly from melons^^. 

8 Bubble equations 

The critical behavior of matrix models can be systematically studied by means of the loop 
equations. In the sequel, we shall generalize the loop equations to bubble equations in arbitrary 
dimensions. Although general bubble equations can be written, we shall present in detail only 
those equations relevant for the leading order in the expansion. 

The bubble equations are a set of constraints satisfied by the partition function of the tensor 
models, essentially, they are the quantum equations of motion. They are derived starting from 
the Schwinger-Dyson equations. The constraints satisfy a Lie algebra. As this Lie algebra is 
somewhat involved, we shall start by detailing it. Only subsequently, shall we relate it to the 
tensor models. 

D-ary tree Lie algebra: While at first sight, this subsection might seem rather technical, 
the ultimate aim is rather simple. Our goal is to construct a composition operator for rooted 
colored {D + l)-ary trees, that is, given two such trees, we can coherently identify another tree 
as their composition. This operation of composition must satisfy certain properties, the most 
important of which is that one can define a set of operators indexed by these trees, which form 
a Lie algebra such that the Lie bracket is compatible with this composition rule. 

From one to an infinity of coupling constants: In this part, we explicitly integrate out D 
of the D + 1 tensors in the path integral to get an effective action and path integral for the 
single remaining color. This effective action contains many effective interaction terms, indexed 
by the possible L'-bubbles of that remaining color. Each such interaction has a (AA)-dependent 
coupling constant. We relax this dependence, rendering the couplings independent. 

Schwinger-Dyson equations: We explicitly describe how the large N limit of the i.i.d. proba- 
bility measure provides a representation of this algebra. 

8.1 A Lie algebra indexed by colored, rooted, D-ary trees 

We have already seen that the leading order graphs are in close correspondence with colored 
rooted D-ary trees. We shall see that the bubble equations we derive are indexed by such trees. 
We first present a set of definitions and properties of the trees. Although technically involved, 
this section is required in order to derive the algebra of constraints satisfied by the partition 
function at leading order (in the large limit). 

Definition 8.1. A colored rooted D-ary tree T with \T\ vertices is a tree with the following 
properties: 

• it has a D-valent root vertex, denoted ( ); 

• it has |T| — 1 vertices that are {D + l)-valent (i.e. each of them has D descendants); 

• it has {D — 1)|7~| + 1 vertices that are 1-valent (i.e. with no descendants). These are the 
leaves we mentioned earlier; 

■^^Fisher's scaling relation occurs in the analysis of critical phenomena of films of finite thickness [65]. In its 
original context, it is a relation among the spin susceptibility exponent, the critical exponent of the spin-spin 
correlation length and the anomalous scaling dimension of the spin-spin correlation function. 
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• all lines have a color index, 0, 1, . . . , -D — 1, such that the D direct descendants (leaves or 
vertices) of a given vertex (or of the root) are connected by lines with different colors. 

We shall subsequently ignore the leaves of the tree, as they can automatically be added once 
the vertices and lines of the colored tree are known. 

Canonical labeling: A colored rooted D-avy tree admits a canonical labeling of its vertices. 
Namely, every vertex can be labeled by the list of colors V = {ii . . . in) of the lines in the unique 
path connecting V to the root ( ). The first color, ii, is the color of the line in the path ending 
on the root (and i„ is the color of the line ending on V). For instance the vertex (0) is the 
descendant connected to the root ( ) by the line of color 0, and the vertex (01) is the descendant 
of the vertex (0) connected to it by a line of color 1 (see Fig. 20 for an example of a canonically 
labeled 3-ary tree with \T\ = 7 vertices). 




Figure 20. A colored rooted D-aiy tree. We have included the leaves in this case. 

Some notation: In the sequel, 'tree' will always mean a colored rooted D-ary tree. A tree 
is completely identified by its canonically labeled vertices, hence it is a set T ={(),... }. We 
shall use the shorthand notation V for the list of labels identifying a vertex. 

Moreover, we denote by i, . . . , i a list of n identical labels i. 

n 

Successor functions: Consider a tree T and one of its vertices: 
V = I , 



l;, . . . , lu, .y, . . . , t.. 



with k ^ i. The colors I, . . . ,k might be absent. The successor of color j of V, denoted 
is the vertex: 



sU{l, ...,k,i,...,i)] 



{I, . . . ,k,i, . . . , i,j) if it exists , 

n 

(/,... ... ,i) ifj/i, 

n 

(/,..., A;) ifj = i. 



if not < 



The colored successor functions are cyclical, namely if a vertex does not have a descendant in 
the tree of color j, then its 'successor of color j' is the first vertex one encounters, when going 
form V to the root, whose label does not end by j. For the example of Fig. 20 we have 



s'[{ )] = 


(0), 


s'[{ )] = 


0, 


s'[{ )] = 


(2), 


s'lm = 


0, 


.^(0)] = 


(01), 


s'lm = 


(0), 


.0[(01)] 


= (01), 




= (0), 


s'm)] - 


= (01) 


Ai2)] = 


(20), 


s'm] = 


(2), 


s'm = 


(22), 
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.°[(20)] = (2), .^(20)] = (20), .^[(20)] = (20), 

s\22)] = (22), s\i22)] = (22), s'[{22)] = (222), 

s°[(222)] = (222), 5^(222)] = (222), s'^[{222)] = {). 

Maximal vertices: We call V the maximal vertex of color i in a tree T if 

V = such that ^^[(^j^^^)] = ( )• 

n n 

In Fig. 20 the maximal vertex of color 2 is (222), the maximal vertex of color is (0) and the 
maximal vertex of color 1 is the root ( ) itself. 

Branches and branch complements: We define the branch of color i ofT, denoted T"*, the 
tree 

r = {iX)\iiX)eT}. 

The branch 7~* can be empty. The root of the branch T*, ( ) G corresponds to the vertex 
(i) G T. The rest of the vertices of T (that is, the root ( ) and all vertices of the form 
{k, U) £ T, k ^ i) also form a canonically labeled tree. Denoted T*, it is the complement in T 
of the branch T^. 

In Fig. 20, the branch of color 2 is the tree 7^ = {( ), (2), (22), (0)}, as all the vertices (2), 
(22), (222) and (20) belong to T. Its complement is 7^ = {( ), (0), (01)}. 

Joining trees: Two colored rooted D-ary trees T and 7i can be joined (or glued) at a vertex 
V £ T. For all colors i, denote the maximal vertices 71 of color i: 

(v^), sVJ(v^)] = (). 

The glued tree, "Ti, is the tree canonically labeled by: 

f(X) foran(X)er, (X)/(y,...), 

{V,Y) foran(y)Gri, 
{V,i,...,i,Z) for ah (y, i, Z) € r. 

rii+l 



T*vTi = { 



This operation can be seen as cutting all the branches starting at V in T, gluing the tree 71 
at V, and then gluing back the branches at the maximal vertices of the appropriate color in 71. 
The vertices of T \ (V) and 71 \ ( ) map one to one onto the vertices of (T*y 71) \ (V), and 
both (y) e r, ( ) G 71 map to {V) G T *y 71, thus \T ★y 71| = |T| + |71| - 1. An example is 
given in Fig. 21, where the leaves are not drawn. 

Equivalence: For any tree T, with maximal vertex of color i, {i, . . . ,i), the maximal vertex 

rii 

of color i in the branch T* will have one less label {i, . . . ,i). One can glue the tree {( ),(«)} 

rii-l 

on (i, . . . ,i) G T*. All the vertices of are unchanged by this gluing, its only effect being to 

rii-l 

introduce a new vertex, (i, . . . ,i) which becomes the new maximal vertex of color i. Subsequently, 

Ui 

one can glue the complement of the branch i, 7"*, on this new maximal vertex: 

T' = (r*(i,...,i){(),(i)}) .,^)r. 
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Figure 21. Gluing of two trees at a vertex T*(2) 71. 



The two trees T and T' have the same number of vertices \T'\ = \T'^\ + \T'^\ = \T\, and the 
vertices of the initial tree T map one to one on the vertices of the final T' ■ As none of the 
vertices of T"* starts by a label i, it follows that the maximal vertex of color i in T' is (i, . . . , z). 

Thus: 

{w) e r, (PF) / (i, y) o («,..., i, T^) G r'. 

Most importantly, it is straightforward to check that the mapping is consistent with the successor 
functions: 

V,W ,W' ^ T' such that s^V] = W ^ s^, [V'] = W' . 

We will say that the two trees T and T' are equivalent, T ^ T' , and extend it by the transitivity 
of ~ to an equivalence relation among rooted trees. An example is presented in Fig. 22. 

The equivalence class of a tree [T] has at most |T| members, all of which are obtained by 
choosing a vertex V = i,j,k,...,l in T and performing the elementary operation ~ on the 
colors i followed by j followed by k and so on up to L 

Properties of colored rooted D-ary trees: We list a number of lemmas concerning the gluing of 
trees, All these properties can be readily understood in terms of the graphical representation 
of the trees. In the sequel, we shall deal with three trees T, 71 and 72. We denote by (z, . . . , i) 

rii 

the maximal vertex of color i in 71 and by {i, . . . ,i) the maximal vertex of color i in 72- The 

proofs of these lemmas are straightforward and all rely on carefully tracking the canonically 
labeled trees under various joinings. They can be found in [81]. 

Lemma 8.1. // {V) = {k, U) e T then {T*v Tif = *u 71, and, for i ^ k, (T *y Tif = V. 
Lemma 8.2. For any three trees, {T *v 71) *v T2 = T *v {Ti *(^) Tz)- 
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T 




Figure 22. Two equivalent trees T ^ V, with V = (j"^ ★(22) {( ), (2)}) +(222) t'^ 

Lemma 8.3. For any two distinct vertices V ^ W £ T, we denote W' the image of W in the 
tree T-kyTi and V the image of V in T*w Ti- Then (T^y 71) *VK' Ti = {T *w T2) *v' Ti- 

Lemma 8.4. Let V £ T and W £ Ti then {T*v Ti) -^w' "Ti = T^y (71 Ti), where again W 
denotes the image of W in the tree T -kyTi- 

Lemma 8.5. We have (7i *( ) Tif ~ (71 ★() Tzf . 

The Lie algebra indexed by colored rooted D-ary trees: We now define a Lie algebra of op- 
erators indexed by the trees. We associate to every tree a variable t-j-, and we denote the 
valency of the root of 71- Consider the operators Cji defined as: 



d 



T 



VeT 



C 



{()} 



N 



D 



T VeT 



d 
dt^' 



where is some parameter (destined to become the large parameter of the tensor model). 
We will not consider the most general domain of these operators. Instead of defining them 
for arbitrary functions /(^t)) we restrict their domain to class functions f{t^j^), with tj-^ = 

Theorem 8.1. When restricted to class functions, the operators Cj- form a Lie algebra with 
commutator: 

[^r2,^Ti]f{t[T]) = ^ ^T2*vTj{t[r\) - X] ^Ti*vT2f{HT])- 
V&T2 VeTi 

The proof of this statement is given in [81] . As it is a straightforward, albeit long computation 
(using all the Lemmas 8.1, 8.2, 8.3, 8.4, 8.5), we do not reproduce it here. 



8.2 From one to an infinity of coupling constants 

In this section, we generalize the i.i.d. colored tensor model with one coupling to an i.i.d. model 
with an infinity of couplings and derive the SDEs of the general model. Our modus operandi 
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here is to integrate all colors save one and then to 'free' the couplings of the operators in the 
effective action for the last color. To get down to business, on integrating all colors save one, 
the partition function becomes: 

Sym(^^) 

where the sum over runs over all connected vacuum graphs with colors 0, . . . , D — 1 (i.e. over 
all the possible L'-bubbles with colors 0, . . . , D — 1) and p vertices. The operators Tr^^ [0^, (j)^] 
in the effective action for the last color are the tensor network operators corresponding to the 
boundary graphs of bubbles, hence provide a natural set of observables of the model. 

Every vertex of is decorated by a tensor cp?^ or , and the tensor indices are contracted 

as dictated by the graph B^ . We denote by v, v the positive (resp. negative) vertices of B^ and 
by lljj the edges (of color i) connecting the positive vertex v with the negative vertex v. The 
interaction operators are written as: 



(d-1 \ 



-D 



i.2) 



where all indices n are summed. Note that, as all vertices in the bubble belong to an unique 
line of a given color, all the indices of the tensors are paired. The scaling in N of an operator is 
determined by the degree of the associated L'-bubble graph: 



2 



where we used equation (4.4) in the first equality. Thus, the effective action for the last color 
(dropping the index D) can be written as: 

_ p-i 

Attributing to each operator its own coupling constant and rescaling the field to T = N ^ (p: 
we obtain the partition function of a colored tensor model with generic potential: 

S{f,T) = ^ t^T^ + ^ts (^"(^)TrB[r,r]. (8.4) 

B 

It is worth noting that, although in the end we deal with an unique tensor T, the colors 
are crucial to the definition of the tensor network operators in the effective action. The initial 
interaction term of the tensor model described a L'-simplex. These later tensor network operators 
describe (colored) polytopes in D dimensions obtained by gluing simplices along all save one of 
their faces around a point (dual to the bubble B). This is in strict parallel with matrix models, 
where higher degree interactions represent polygons obtained by gluing triangles around a vertex. 

When evaluating amplitudes of graphs, obtained by integrating the last tensor T, the ten- 
sor network operators act as effective vertices (for instance each comes with its own coupling 
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constant). It is still more convenient, however, to represent the Feynman graph of the path 
integral (8.4) using graphs with D + 1 colors. The effective vertices are the subgraphs with 
colors 0, . . . , D — 1, and encode the connectivity of the tensor network operators. 

The partition function (8.4) provides a natural set of observables for the model, the multi- 
bubble correlations, defined as: 

TVh(„ [f, T] Tr6(„ [f, T] • • • TVb(,) [f, T]) = fl C^^^" (^"(^«^] Z. 

When introducing an infinity of coupling constants, we did not change the scaling with 
of the operators. The graphs G contributing to the connected multi-bubble correlations are 
connected vacuum graphs with D + 1 colors and with p marked subgraphs corresponding to the 
insertions Trg^^^^ [T, T]. Taking into account the scaling of the insertions, the global scaling of 
such graphs is: 



TrB(„ [T, T] TrH(„ [T, T] • • • TVb^^JT, T]^ 

where we use Lemma 4.6. 

In the large limit, the connected correlations receiving contributions from graphs of de- 
gree (melonic graphs) dominate the multi-bubble correlations. All their bubbles are necessarily 
melonic, in particular, the insertions Trg^. JT',T]. As we have seen in Section 7.1, the melonic 
D-bubbles (i.e. melonic graphs with D colors) are in one-to-one correspondence with colored 
rooted D-ary trees T- The tensor network operators (8.2) of melonic bubbles can be written 
directly in terms of T. When building a D-bubble starting from T, each time we insert a melon 
corresponding to a vertex V € T we bring a T and a T tensor for the two external points of the 
melon. We denote the indices of T by ny and the ones of T by ny- We get: 

/ D~i \ 



Tvb[T,T] ^ Tvr[f,T] = n T^v% n '^^ 



where Sj- is exactly the colored successor function defined in Section 8.1. As this operator 
depends exclusively of the successor functions, it is an invariant for an equivalence class of trees: 
T~ T' Tr7-[r, T] = TiC'j-'[T,T]. Hence, the action and the partition function depend only on 
the class variables tjT-j = ^7-/^7- ^t- Taking into account that the melonic bubbles have degree 
(and redefining the coupling of the tree T = {( )}), the action may be written as: 

S{f, T) = Y,tT Trr [r, T] + S^{f, T) , 

r 

where S"" corresponds to non-melonic bubbles. 
8.3 Schwinger— Dyson equations 

Consider a melonic bubble corresponding to the tree 7i with root ( )i. We denote: 

'^Sn = n n 



The SDEs are deduced starting from the trivial equality: 



E 



6 




0, 
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which expands as: 



E 



' n ^...] ^ 



V2^( )i 



n 



,yieri\{( )i,V2} 



T- I 8^^- 



^Vieri\( )i 



(8.5) 



The second line in (8.5) represents graphs in which a line of color D on a melonic bubble connects 
the T on the root melon ( ) to a T on a distinct melon V2. Hence, it can not be a melon (see the 
end of Section 7.1). The last term represents a melonic bubble connected via a line to a non- 
melonic bubble (coming from ^§jr). Thus, it cannot be a melon either. Taking into account 
that we have one line explicit in both graphs, (hence a factor A^^*^^^^^), and that the scaling of 
Tr7-[T, T] is N^~^, we find that in both cases the correlations scale at most like 

1 



.(^...'j (D-l)!. 



(8.6) 



The first term in (8.5) factorizes over the branches 7i of 71. We denote by ( )i^j the root of 
the branch Ti- For a non-empty branch T*, recall that the vertex sij-[{ )i] = (i) G T maps to 
the root ( )i^j G T*. For each branch, we evaluate: 



. ^ Oi Oi Oi sL- [()i 



)l 



N 



= Sj^i if not. 

4^ [Oil l^tj-i]-^li)i.i] Ci.i 

Thus, denoting the valency of the root ( )i G 7i by we get: 

D-l 



i=0 



The third term in (8.5) simplifies: 



D-l 



D~l 



\" 8- 5^ = \" 8- T\ 8 T\ 8 ^7-i\()i.r\v 

. ^ ^ ^ .-—O ' Ti' "■•'1' n 



=0 



D-l 

^n^ . , n' "?i,n "n,n ^n,n 

Hence, the SDEs, for every rooted tree 71 with non-empty branches, may be written starting 
from the root as: 

nD-\r.\ / Y[ Trr.[f,T]) - N^-^Y.'r E (Trr.,ri[T,T]) = (...>, 

\ i=0 
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where ( • • • ) denotes the non-melonic terms of (8.6). Taking into account the definition of £7-^ 
in (8.1), we spot that: 

Lr,Z = {■■■) hm (N-^'^Cr.z) =0, VTi. 

W— >oo \ Zj J 

Recall that Z = e^^^^*^^ depends only on class variables t^j^. At leading order in only 
melonic graphs contribute to the free energy FUb), hence lim F{tB) = Foo{\i\T\\)- The SDEs 
at leading order imply: 



i=0 



The most useful way to employ the SDEs is the following. Consider a class function Z = 
g-Af F satisfying the constraints at all orders in A^, C-j-jZ = 0. Its free energy in the large N 
limit is Foo{tm) = lim F{t!'J-^) and respects (8.7). As a result, Foo(im) = -^oo(im)j that is, 

the A — )• cx) limit of Z and Z coincide. 

Note that an SDE at all orders can be derived for the trivial insertion: 



E 



_6_ 



p 



where p denotes, as usual, half the number of vertices of the bubble B. The above operator, 
which at leading order is C(^-j, should be identified with the generator of dilations [55]. 



9 Non-trivial classical solutions 

It seems unlikely that we shall ever have at our disposal an exhaustive catalog of classical 
solutions to these models, since most quantum field theories are in a similar position. Quite the 
opposite situation is more prevalent; it is non-trivial, especially for scalar theories, to have any 
solutions at all in closed form (other than the trivial solution). 

Here, we shall formulate an ansatz to extract classical solutions to colored tensor models. On 
the one hand, the ansatz is rather general in the sense that it appears to be insensitive to the 
specific details of the covariance and vertex kernels, although we shall only apply it explicitly 
to the i.i.d. model. On the other hand, it is just an ansatz and so for a specific model it only 
captures a small fraction of the solution space. It is based upon a splitting of the tensor model 
structure into a background and a foreground. The PJ-factorizations of Section 4.7 enable this 
mechanism. 

In fact, we are not just interested in acquiring non-trivial classical solutions. In the sequel, 
once we have found explicit classical solutions for the background, we shall be interested in 
analyzing the quantum theory resulting from fluctuations around this classical point. It emerges 
that a fruitful area of research lies not just in analyzing the most generic fluctuations but also 
more restricted scenarios. 
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9.1 Ansatz set-up 

The ansatz has three distinct steps: splitting the degrees of freedom, solution finding and 
quantum analysis. 

Splitting degrees of freedom: Sticking with the i.i.d. probability measure, as usual, the first 
step we take is to separate the degrees of freedom of the tensor model into two types, background 
and foreground. We recall from Section 6 that a neat way to do this splitting is to utilize a PJ- 
factorization. This packages the degrees of freedom of the tensor model into sets pertaining to 
jackets and patches. From this, one can split the jackets (and patch) into background jackets 
and foreground jackets. Let us take D = 3 as an example. Every PJ-decomposition contains 
one jacket and one patch. Thus, one has two choices, once a given factorization is chosen: either 
the jacket forms the background and the patch the foreground or vice versa. 

Classical solution: As a next step, one attempts to find solutions to the theory in the sector 
where one has trivialized the foreground. In essence, one looks for solutions to a lower rank 
tensor theory albeit still with D + 1 colors. In realistic terms, of course, we are not ready to 
examine all such sectors. Our algebraic techniques are suitably advanced only in the areas of 
vector and matrix models, that is, backgrounds consisting of a patch or a jacket. More involved 
backgrounds are still out of our reach. For L> = 3 or 4, however, these simpler backgrounds are 
sufficient to exhaust the ansatz given above. 

Effective quantum theory: The power of the splitting becomes amplified once one returns to 
the quantum theory. It allows one to parse systematically through various types of quantum 
fluctuation. 

The first and most general fluctuation corresponds simply to expanding the original tensor 
about our new solution. In this way, we do not depart from our original theory, we rather view 
it in a new and hopefully revealing light. 

It is also interesting to consider less general quantum fluctuations. Perhaps the least gen- 
eral corresponds to restricting them to consist of only foreground degrees of freedom. For the 
i.i.d. model, as we shall see shortly, these quantum fluctuations decouple from the background. 

There exist also some intermediate scenarios. We can choose to allow a mixture of background 
and foreground fluctuations. This allows for some back-reaction onto the background. In fact, 
once one has a specific background, one can analyze fluctuations corresponding to an arbitrary 
collection of jackets. These jackets need not even be in the original PJ-factorization. Such 
a jacket contains degrees of freedom which are partially foreground and partially background 
with respect to the original slicing. 

Remark 9.1. We refer the reader to [63] and [105], where a systematic analysis of the quantum 
fluctuations around classical solutions of the Boulatov-Ooguri model in 3 dimensions was con- 
ducted. In both cases, the background is a patch, while the foreground is a jacket. Due to the 
non-trivial covariance of the model, one finds that even purely foreground fluctuations do not de- 
couple from the background. Interestingly, these jacket fluctuations can be re-interpreted as the 
quantum degrees of freedom of a non-commutative massive scalar field theory [69]. Moreover, 
expanding the action to quadratic order in the fluctuating fields, these jacket fluctuations along 
with pure patch fluctuations provide a basis for the full quantum fluctuations. In other words, 
at quadratic order, the quantum fluctuations are diagonalized in terms of a PJ-factorization. In 
the full action of course, interference occurs due to higher order interaction terms. 

Remark 9.2. These classical solutions are often referred to as instantonic solutions in the group 
field theory literature, although this is a slight misnomer. Genuine instantons in field theory lie 
at minima of the potential, meaning that the value of the action at these solutions is less than 
that for the null configuration. Thus, they contribute heavily to the path integral, although 
this would be missed at low orders in the naive perturbative expansion. It makes sense then 
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to treat the quantum theory as a perturbation about these instantonic solutions, which provide 
the dominant contribution to the semi-classical physics. We are not in this scenario here. Our 
classical solutions are not minima of the potential. To see this, one can appeal to the derivation 
of the effective action in Section 8.2. One is left with a single color, say, once all the other colors 
{1, . . . ,D} have been integrated out. As we saw in (8.3), all the coupling constants (AA)-'/j! 
are positive. Thus, the true minimum of the effective potential is cp^ = 0. This is also the case 
for the simplest matrix models, but one may alter the coupling constants in the effective action 
to get to a different theory with true instantonic contributions. In other words, some of the 
coefficients need to become negative. To see that we are not hiding anything by integrating out 
these fields, note that the positivity of the coefficients emerges even after integrating out just 
one color, D say. The potential is once again positive definite and the minimum is at: = 
for ie{0,...,D-l}. 



9.2 Patch background 

Let us consider the classical equations of motion of the {D + l)-colored i.i.d. model. There are 
2D + 2 such equations: 
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which expand as: 



D(D-l) ii 



and their conjugates. In its present incarnation, we have not lost the manifest division of 
degrees of freedom among the jackets of a PJ-factorization as in (6.1). Thus, having decided on 
a background, we can pick solutions which trivialize the foreground sector. The simplest case 
is where the background corresponds to the patch in a PJ-factorization. To possess a patch, D 
must be odd. We choose the specific PJ-factorization constructed earlier in Proposition 4.7 and 
consider solutions of the form: 



$1 



where 5nj is shorthand for a product of identity matrices associated to the background jackets 



in our PJ-factorization. This is our ansatz. These functions must satisfy conditions imposed by 
the equations of motion: 



A 



_D{D-1) 



N- 



(D-l)/2 

n ■ 

j=0 



D+l 



(9.1) 



where we have introduced vectors /* = /V, d+i ■ One can explicitly realize functions satisfying 

these equations but we shall leave our treatment in a rather generic state here. Now, one can 
expand the field theory about this non-trivial classical solution as = + , leading to: 
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is the action with respect to this background ^> and we have used the equations of motion (9.1) 
to remove the terms linear in and (p^ . Most clearly, we see that we now have a field theory 
with a non-trivial covariance C$ dependent upon the background $. This could be of use when 
one wishes to devise a renormalization prescription for the theory [105, 116]. Of course, one 
can restrict to purely foreground fluctuations. However, in this general form one does not see 
anything particularly striking. We refer the reader once again to [63, 105] for specific results 
relating such quantum fluctuations to matter field theories and the Hamiltonian constraint 
operator of quantum gravity in three dimensions. 

Remark 9.3. As already mentioned, one need not take the background to be a patch. The next 
simplest case, corresponds to a background consisting of a single jacket, that is a single Riemann 
surface, and is possible in any dimension D. One needs to solve a set of coupled algebraic matrix 
equations in that case. 

10 Extended discussion and conclusion 

We have by no means attempted to compile a historical or complete review of higher-dimensional 
tensor models. On the contrary, we wished to commence from what we believe are the core 
mathematical precepts upon which the whole theory is based. Thence, we felt it was perhaps 
more beneficial that the review developed in a rather linear deductive style, obviously guided by 
particular prescribed motivations. While these of course are bound to be somewhat idiosyncratic 
but they are by no means uncommon; they are some of the basic questions one should expect 
a quantum field theory to answer. With these provisos, we have pushed right up to the cutting 
edge. It does mean, however, that we arrive at a rather singular focus, perhaps at the expense of 
other topics. With this in mind, we present here a short summary of results that lie at a slight 
tangent to, although by no means disconnected from, the main text. We divide this section in 
three parts. The first deals with further developments on critical behavior(s) in colored tensor 
models. The second details the connection between colored tensor models and spin foams. The 
third and final part draws the conclusions of our work. 

10.1 Critical behavior in more complex models 

Now that the main ideas behind the large N behavior of colored tensor models are in place, it 
is only natural to apply them to particular models describing specific physical situations. In 
particular a very active research direction is to identify more complex critical behavior in various 
models. 

The Ising model on a random lattice in higher dimensions. One of the most impressive early 
successes of random matrix models was the analytic solution by Kazakov et al. [94], followed 
by [39, 41], of the Ising model on a random planar surface. In particular, the solution provides 
analytic proof that the system undergoes a third order phase transition between the low tempe- 
rature magnetized phase and the high temperature non-magnetized phase. A similar study [34] 
can now be undertaken in arbitrary dimension for the Ising model on a random melonic lattice. 
As we explain below, in the continuum limit, the spin system does not exhibit a phase transition 
at finite temperature, which is in agreement with numerical investigations. 

The two-matrix model [94] used to describe an Ising spin system generalizes to a two tensor 
model, both colored, (X*,y*) (and (X*,y*)) defined by the partition function: 
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D 

s{x, X, y, y) = y y [xl xl + yi yi - cxi yi - m xi 



2=0 U. 



D _ D 

X 



fc j=0 ^ 1=0 

The Feynman graphs of the action (10.1) are identical to the ones of the i.i.d. model, up 
to the fact that a vertex now involves either X or y tensors (both coming with positive and 
negative orientations). Tensor indices kij and colors are preserved along the strands, and the 
propagator between the X and Y sectors is 



C 



1 c 



1-C2 \C 1/ ■ 



Lines in a graph can join a field X (resp. Y) to a field X (resp. Y), with weight 1/(1 — (?), or 
a field X (resp. Y) to Y (resp. X), with weight c/(l — c^). This is the only coupling between 
X-tensors and y-tensors. As in the i.i.d. model, the free energy organizes in powers of 1/A^, 
with melonic graphs dominating the large limit, so that 



s{g) xP-x^-yPyyPy—-^+0{N-^), 



melons Q 

with px (resp. Px, Py and Py) the number of X (resp. X Y and Y) vertices, L is the total number 
of lines, Lj^y (resp. L^y) is the number of lines from a X to a y (resp. from a X to a y) 
and s{Q) is a symmetry factor. We have L = {D + l){px + Py)-, Px + Py = Px + Py- 

The X-vertices and y- vertices represent the Ising spins on random lattices. One, say X, 
mimics the random world of the up spins, while the other the random world of the down 
spins. On a graph Q, the total number of vertices n is the number of Ising spins. Through the 
propagator, two parallel spins are coupled with the same weight 1/(1 — c^) and two anti parallel 
spins with c/(l — c^), from which one can extract the Ising coupling /3J where /? is the inverse 
temperature. Taking c E [0, 1] corresponds to the ferromagnetic model J > 0. Setting x = x 
and y = y, the free energy of the tensor model equates the grand canonical partition function 
of the Ising model on random lattices: 

Z{^l, (3 J, Ph) = Y e""^^n(/3 J, m, 

n 

Zni/3J,l3h)= Yl s(g) ^ e'^^^<'.^)^'^^+'^'^^»^\ 

g,2{Px+Py)=n Si=±l 

Zn is the canonical partition function for the Ising system on random lattices with a fixed number 
of spins n, h is the magnetic field and fi the chemical potential. The relationship between the 
Ising parameters and those of the tensor model is: 

/ 1 \i/2 /t\^/^ 



e = c, 



As the model depends on two complex tensors, there are four types of 2-point functions, denoted 
^XX), ^yy), (Xy) and (yx), with which one must deal. In the melonic sector, the self- 
energy factorizes into convolutions of 2-point functions. Some algebra leads one to a system 
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of equations for the two point functions (at zero magnetic field), {XX'^ = (YY^ = U and 
{XY) = {YX) = V: 

U - cV - g{U^^^ + V'^+^) =1, V - cU - gUV{U^-^ + V^~^) =Q, (10.2) 
The solution for the noninteracting model (g=0) at arbitrary temperature is: 

c/(c,o) = -^, y(c,o) 



which is a boundary condition for the system of equations (10.2). For arbitrary dimensions, one 
can obtain a parametrized solution of (10.2). To see this, we introduce the following pair of 
variables: 



1 

Uil-gU^y 



and rewrite equations (10.2) as: 

c = (1 - z'')w - z^'w'', gi/^ = z{l - z^)(l - w^). 

For sufficiently small values g, the free energy (related to z in the new parameterization) is 
analytic. We now send the lattices to criticality. We obtain a critical curve gdc), on which the 
free energy F{c) = F(c,gc{c)) describes the Ising system in the large volume limit as a function 
of the temperature. 

The critical curve is defined by the vanishing of [dg / dz^ , so that the expansion of g in terms 
of z only begins at the second order. This equation has four solutions ^^i, 2, 3,4(^^)1 and one can 
show that the physical solution zi{w) is isolated form the others at finite temperature, hence 
the free energy is analytic for all < c < 1. In the infinite temperature limit c — ?■ 1, the physical 
root zi(w) collapses to and the free energy acquires a non-analytic behavior. The Ising spin 
system on random melonic lattices undergoes a phase transition at infinite temperature. 

Dually weighted models. A second model, which we shall discuss here, consists of a class of 
colored tensor models with a modified propagator which allows us to associate weight factors to 
the faces of the graphs [25]. Such models correspond to dynamical triangulations in three and 
higher dimensions with generalized amplitudes. We present a model which undergoes a third- 
order phase transition in the continuum characterized by a jump in the susceptibility exponent. 
This model is inspired by the dually weighted matrix models introduced by Kazakov et al. [26, 
95, 96, 97, 127]. The dually weighted colored tensor model in dimension D is defined by the 
partition function: 



A ^ . A ^ -. 

j\rD(D-l)/4 X/ n jY^(D-l)/4 X] n "^iti ■ 

n i=0 n i=0 

Using the counting of faces of a colored graph established for the i.i.d. model, the amplitude of 
a graph can be written: 

~,T\P/ 
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where 2p|p) is the number of vertices of the p'th face with colors (ij) and is the matrix 
transpose C. The fuh connected 2-point function of the dually weighted model: 

is independent of the color i and factorized along strands. Surprising as it may seem, the self- 
consistency equations in the melonic sector can be solved analytically for P [25]. A particular 
model one can consider is defined by a covariance C such that: 

TT[{C^Cy] =Nq-^. 

Such a choice corresponds to the DT amplitude: 

A{g) = e''D-2Nn-2->^DNu -Q q-P^ (10.3) 



where the product is over all {D — 2)-dimensional simplices (bones) of the triangulation, with qi 
being the number of L'-simplices to which the bone i belongs. The DT amplitude (10.3) was 
studied via numerical simulations in [30], and more recently in [102]. 

The self-consistency equation (and its physical initial condition) appear in the melonic sector 
(denoting gdgE = U) as: 



U = S{^, z{9, U)) , 5(/3, .) = g ('J) 



JJ^ D(D+\) 2 

<9,U) = 5(/3,0) = 0. (10.4) 

We shall denote the solution of the equation (10.4) by g{j3, U). The series 5(/3, z) has radius of 
convergence: Zb = ^ for all values of /3. 




U 



Figure 23. A schematic representation of the phase transition and the low /? and high f3 regimes. 

The self-consistency equation (10.4) admits two regimes corresponding to two distinct phases 
in the continuum limit. Either, for small /3, the curve g{f3, U) achieves its maximum before 
exiting the analyticity domain of S, in which case, one recovers the usual critical behavior: 

^non-analytic ^ {d ~ 9c)^ ■ 

Alternatively, for large enough /3, the curve g{(3, U) exits the analyticity domain of S before 
reaching its maximum, in which case, we recover a non-analytic behavior due to S{(3, z) itself: 

-E'non-analytic ~ {§ — QhY^ ^ ■ 
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The critical value /3c, where the phase transition occurs can be determined numerically: Pc ^ 
1.162 for D = 3, l3c^ 1.216 for D = A, and Pc » 1-134 for D = oo. A schematic portrait of the 
transition is presented in Fig. 23. The critical behavior is different at /3 = f3c, namely [25]: 



-E'non-analytic ~ " 5c) ''"^^ . 

In conclusion the susceptibility exponent, defined by -Bnon-anaiytic ~ (5 — gcf'~'^ ^ is: 
\ for ^ < /3c, 

7 = \ T = /^c, 

Pc+ 2 

\-p for/3>/3c. 

We observe that for a large range of values of /3, that is, for /3 < /3c, universality holds: the 
critical exponent 7 is independent of /3. On the contrary, for /3 > /3c, one has a one-parameter 
family of different critical behaviors. A possible interpretation of such unusual behavior is that 
for /3 sufficiently large the new term starts behaving as a non-local or long-range interaction, 
for which universality is not expected to hold. This /3-dependent behavior, which we encounter 
here, is reminiscent of certain models of branched polymers [4, 29]. 

One can further show [25] that the phase transition is third order. This result represents 
the first time that a phase transition in the continuum limit is accessed by analytical means in 
dimensions higher than two. Furthermore, it is the first time that a third-order phase transition 
is observed in DT, which opens up the possibility of obtaining a continuum limit with 2 < 
dn < 00, without the causality condition employed in causal DT [8, 9, 10, 12]. The link to 
the numerical results of [102] deserves further exploration, either by pushing the simulations to 
larger k_d-2 and larger volumes, or by trying to extend the analytical tools to finite hd-2- 

There are of course many interesting open questions regarding the critical behavior of colored 
tensor models. One would like for instance to find multi-critical points generalizing the classical 
multi-critical points of matrix models [53]. Also a better understanding of the continuum limit 
can be gained via some sort of double scaling limit including contributions from sub-leading 
terms in the 1/A^ expansion (corresponding to a finite kd-2 in the DT interpretation). Third, 
in two-dimensional matrix models the KPZ correspondence [47, 57, 98] relates statistical systems 
on random and fixed lattices and led to early computations of novel critical exponents [60] . One 
would like to find the appropriate generalization of this correspondence in higher dimensions. 

Renormalization and colored tensor models. A direction of research in which the 1 /N expan- 
sion of the colored tensor models is very consequential is the renormalization of non-identically 
distributed tensor models. We shall briefly draw the lessons one learns from the results presented 
in this review. When studying non-i.i.d. models, the divergences regularized by the re-scaling 
with of the coupling constant become UV divergences in the amplitudes. Such divergences 
must be dealt with using renormalization group tools. A similar situation appears in matrix 
models where the renormalization of non-i.i.d. models [78, 79, 88, 89, 117] is intimately related 
to the 1/N expansion. As some non-i.i.d. tensor models are UV complete [22, 56, 80], one hopes 
that a non-i.i.d. tensor model can provide a non-trivial complete quantum field theory. 

Renormalization is a recipe relying on three intertwined ingredients [116]: a scale decom- 
position, a power counting, and a locality principle. In the case of matrix models, both the 
locality principle and the power counting can be read off already for the i.i.d. models. In the 
1/N expansion: 

• planar graphs with one external face dominate. This is a key property, since in non- 
i.i.d. models, only planar graphs with one external face can be subtracted by adding local 
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counter-terms in the action (the original trace interaction of the action are planar with 
only one external face). The 'locality' of usual quantum field theories becomes, in matrix 
models, 'planarity with one external face'. 

• under the appropriate scaling of the coupling constant, the planar graphs with one external 
face are uniformly divergent with N yielding an infinite family of uniformly divergent 
graphs. In non-i.i.d. models, this leads to a non-trivial subtraction and running of the 
parameters. 

In non-i.i.d. models the propagator must scale for large matrix indices exactly like the ap- 
propriate power of in order to conserve these two features. Planar graphs will dominate 
any model, but any other scaling of the propagator leads either to non-renormalizable or super- 
renormalizable models. 

The melonic graphs act in higher dimensions as the planar graphs of matrix models. They 
dominate the 1/N series having the largest number of faces for a fixed number of vertices. As 
such they will dominate any model. In order to write a sensible non-identically distributed model 
one must propose a propagator with the appropriate scaling for large indices, which render the 
melonic family uniformly divergent with respect to the cutoff. Such divergences can then be 
subtracted by adding melonic counter-terms in the action. 

One can ask if it is possible in any way to subtract the melonic family. That is, is it possible 
to write some tensor models dominated by non-melonic configurations? It is premature to give 
a definite answer to this question, but we doubt that this is possible. Non-trivial topologies are 
generated by (core) graphs with very few vertices. For instance in D = 3, already at 6 and 8 
vertices one encounters not only the MP^ manifold but also pseudo-manifolds. If one subtracts 
the melons it is very likely that non-trivial (non-manifold) topologies will dominate the effective 
behavior of the series. 

10.2 Tensor models and spin foams 

Each graph of a colored tensor model represents a triangulation. If the base group is SU{2), 
for instance, the faces of the graphs are then labeled by spins (which one sums to recover the 
amplitudes of the graph). Such a triangulation with spins associated to its {D — 2)-simplices is 
called a 'spin foam'. The spin foam models are then particular realizations of tensor models on 
fixed graphs. 

A brief prelude - the colorless world: When concentrating on particular graphs, one can 
of course chose a preferred graph on which to conduct a study. As one need not deal with 
the full series generated by a path integral, in this setting, one can consider also non-colored 
graphs. In this case, the connection to cellular topological spaces is by no means direct, since 
the cellular structure is no longer intrinsically encoded in the graph. Actually, these graphs 
are only 2-complexes, possessing just 0-, 1- and 2-bubbles. In other words, the higher simplices 
alone have a matching prescription. This leaves a lot of implicit freedom in how one identifies 
lower-dimensional simplices as one moves to ever higher dimensions. Essentially, there exists 
a class of diagrams (dominant both by weight and entropy) which cannot be immediately related 
to simplicial pseudo- manifolds, rather only after sufficient barycentric subdivision [50]. This 
barycentric subdivision can not be encoded in the Feynman rules of a tensor model. 

As a result, one is not in a position to perform a expansion, extract the dominant order, 
examine quantum equations of motion; it appears that the colors are indispensable for this 
purpose. Despite that, however, uncolored models have a place in the scheme of higher tensor 
models since in a certain sense they allow one to test whether or not certain properties depend 
intrinsically on the stricter topological control one exerts in colored models. For example, the 
Borel summability of some models does not necessarily require that the model is colored [70, 107]. 
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Unfortunately, even this result requires a modification of the theory by adding an extra term 
(the pillow term) which changes the perturbative power-counting. 

The 'projector' models: Undoubtedly, a motivation for the study of tensor models comes from 
a discrete approach to quantum gravity known as spin foam quantization [66, 111, 112]. This 
particular brand of theory relies on a simplicial (pseudo-)manifold, upon which the gravitational 
degrees of freedom are first placed and then quantized. 

Three-dimensional gravity is ideally suited to this procedure, since it may be recast a topo- 
logical gauge theory. In this case, the absence of local degrees of freedom facilitates the passage 
to a discrete theory without any serious anomalies, in the sense that the local gauge algebra is 
preserved at the simplicial level. Quantization in a lattice gauge theory sense is then immediate. 

This case cannot be made for gravity theories in higher dimensions. There, the simplicial 
manifold implies a severe truncation of the local degrees of freedom. Also, there appears to be no 
way to preserve diffeomorphisms. One might hope to still possess some seed which would lead 
to diffeomorphism symmetry in an infinite refinement limit, or one could hope to regain these 
degrees of freedom by summing over infinitely many triangulations. An initial step for the first 
method is to find a coherent refinement strategy while for the second, it is to find a summing 
prescription. Colored tensor models offer a canonical summing prescription (alternative, less 
canonical, summing prescriptions have also been proposed [118]). 

All current spin foam models place weights just on the 0-, 1- and 2-bubbles. In particular, we 
are talking about the Boulatov-Ooguri [38, 110] and EPRL/FK models [21, 61, 62, 68, 106]. To 
give a flavour, let us give the general flavour of these non-i.i.d. models with specific projectors 
provided for the Boulatov-Ooguri (BO) class of models: 



du{(t),(t)) = YldfiT^Bo{^\^')e ^4«)^{"}nf=o'^{«j ^4o)E{S}nf=o^=,^ 
V^o^ =y6^.., aBO(^) ^ {D-l)iD-2) ^ 

rrii 

One can see that the Boulatov-Ooguri projects out a diagonal component, while the value 
of is chosen to make a 1/A^-expansion meaningful. From the tensor model point of view, 
they represent some particular choices of non-i.i.d. covariance, while the vertex kernel remains 
unchanged. Note that in the literature such models are sometimes seen as a modification of 
the vertex while the propagator remains unchanged. Of course, by dressing the field itself with 
the square root of the propagator one can shift the non-triviality from the covariance to the 
vertex kernel and vice versa. The two formulations are perfectly equivalent, but we favor the 
former. The rationale behind this is that for the usual (/>^ model, such a rescaling leads to an 

4 

action with 6(pi — P2) as a propagator and 6(pi + P2 + P3 + P4) FI , ^ interaction kernel. 

a=l ^Vl+vn? 

One can of course establish renormalizability of the c/)^ model in this setting, but it is more 
difficult to prove that the vertex kernel is in fact reproduced under the renormalization group 
flow. In addition and more importantly, since the quadratic part is just the identity, one might 
be tempted to search for the semiclassical behavior in the large (semi-classical) p limit. In the 
case of the (/)^ theory, this is certainly a misguided endeavor, since the renormalization group 
flows in the opposite direction. For spin foams, this argument at the very least cautions against 
prematurely drawing conclusions from results obtained in the large spin limit. 

The Boulatov-Ooguri tensor model assigns a BF theory weight to each graph [31]^^. There- 
fore, in the 3-dimensional case, it is appropriate for describing gravity on a simplicial pseudo- 
manifold. One may wonder why one would care about summing over graphs in this case, given 



^^BF theory is a topological gauge field theory dependent on two fields: a g-valued (D — 2)-form B and a g- 
valued 1-form connection w, which takes part via its curvature F[w\. g is the particular algebra of interest. Its 
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that one has a well-defined quantum theory on a graph-by-graph basis. One may appeal to the 
2-dimensional case for a certain amount of motivation. The Einstein-Hilbert action in that case 
is a topological action; it is sensitive only to the Euler characteristic of the underlying manifold. 
Upon quantization, however, the theory acquires the non-trivial property that quantum fluc- 
tuations may induce tunneling between topologies. Thus, a sum over topologies is appropriate to 
capture the true quantum nature of the theory. Moreover, in the presence of a cosmological term 
the theory becomes highly non-trivial [57] . While such a direct motivation is less persuasive in 3 
dimensions, the sum over topologies can induce non-trivial effects. For the colored version of the 
Boulatov-Ooguri models, one can again construct a topological 1/N expansion and again melonic 
graphs are the sole contributors to the leading order [85, 86, 90]. One can hope for a critical be- 
havior reminiscent of i.i.d. models, but one should keep in mind that cutoff effects, although not 
modifying the power-counting of graphs, might influence the (previously combinatorial) finite 
coefficients of the graphs in the melonic series. Various bounds and amplitude evaluations have 
been established for Boulatov-Ooguri like tensor models [20, 23, 24, 35, 36, 37, 45, 67, 107]. 

The EPRL-FK model assigns a constrained BF amplitude to each graph, intended to mimic 
the gravitational theory for D > A. Here, of course, one hopes to use the full power of tensor 
models to regain a classical gravity on a differential manifold (symmetries and all) in the con- 
tinuum semi-classical limit. A expansion has not yet been constructed for these models 
and constitutes an open research question. Some intriguing results on graphs amplitudes have 
been established in the literature [101, 115] for these models. 

Bahy spin foams: A recent development in the spin foam literature is the advent of so-called 
'baby spin foams' [14]. These are spin foam models utilizing finite groups (for example, Z^v) 
rather than the Lie groups (as used in Boulatov-Ooguri and related models). The hope is that 
the simplified group structure might allow one to perform more complicated tasks within the spin 
foam formalism: to perform calculations with a large number of simplices, at least numerically 
and to utilize methods from statistical lattice theories to coarse-grain spin foams [58]. Moreover, 
these models provide a nice middle ground between tensor models and spin foams, since both 
the i.i.d. and non-i.i.d. potentials detailed earlier generate spin foam amplitudes based on the 
group Zat or a compact Lie group with a cut-off. 

Symmetries and Ward-identities: In analyzing various colored tensor models one needs to 
address the question of the symmetries of the action. Indeed, they are hugely consequential for 
the renormalizability of the theories. One can adopt a more algebraic approach and relate the 
symmetries of the tensor to an underlying n-ary algebra [122, 123, 125] of functions. In this 
setting, every given tensor configuration corresponds to an algebra whose structure constants 
are given by the tensor. In a more quantum field theoretic perspective, the issue of symmetries 
is neatly clarified in the colored models. Indeed, as the tensors have no symmetry properties 
one can perform local shifts in various arguments without affecting the others. This allows 
the analysis of such symmetries at the classical [18] and quantum [19] level. Moreover, the 
colors have proved instrumental in the identification of a peculiar symmetry of projector models 
under shifts of the arguments corresponding to 0-simplices (vertices) of the triangulation [15, 73]. 
Understanding in more depth the symmetries of tensor models is a very active research direction. 

action is of the form: 



In the case of 3d gravity is a smooth manifold, g = su(2) and k — — IGttG. For such a theory, one replaces A4 by 
a simplicial presentation A of the same topology and the smooth fields by counterparts which are distributional 
on relevant sub-simplices. Thereafter, one arrives at an analogous discrete theory, which nonetheless encodes 
the relevant information of the continuum theory, since both capture just topological features of the manifold in 
question. In the example given in equation (10.5), g = Z„, but one may formulate the tensor model for gravity 
inspired gauge groups. 
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There are many other topics in tensor models to which this review does not do justice. From 
alternative models [130] to the intimate relation with non-commutative geometry via a Fourier 
transform [16] (which naturally sheds light on the geometry of the quantum space time [17]); from 
the emergence of matter [54, 63] on a non-commutative background to collective behaviors [113] 
or topological polynomials for tensor graphs [129]; these and may other topics are left for the 
dedicated reader. 

10.3 Recapitulation 

Alas, we have arrived at the moment to close off this review of colored tensor models. The 
preceding sections have been largely focused on adapting and applying methods well-known in 
statistical physics and matrix models. We paid less attention to the physical and mathematical 
applications of the colored tensor model framework. For sure, there is a multitude of toy-model 
scenarios of physical theories, which have been known for quite some time to be captured by 
tensor models, but as for real-world examples .... Of course, we hope that rapid progress in 
this area will rectify this situation sooner rather than later. All the same, we expect that the 
core mathematical techniques, which we presented here, will continue to act as a linchpin for 
subsequent developments, no matter what direction they might take. With this in mind, we set 
down our roots in the theory of probability measures, stripping many of the implicit physical 
connotations one inherits from a specific application. This allowed us to concentrate on the core 
principles without the necessity to detail the physical interpretation at every step. 

Thereafter, we provided a brisk review of some familiar techniques employed in multi-matrix 
models with a view to extending them to higher dimensions. This set the scene for the subset of 
quantum field theories to which this review is devoted, the (-D + l)-colored tensor models. From 
a perturbative point of view, the central feature is that the Feynman graphs encode, as succinctly 
as possible, D-dimensional simplicial pseudo-manifolds. The form of the tensor models permits 
a sharp division of labor: the colors are solely responsible for the cellular topology, while the 
index structure is responsible for how the sub-cells are weighted in each cell complex. Since this 
has so much potential, we embarked on a detailed study of a subset of structures embedded in 
the D-complex. The bubbles obviously play an important role no matter what the choice of 
tensor, covariance or vertex kernel. One needs them to grasp the processes of dipole creation and 
deletion, which constitute graph manipulations, including, as a subset, the homeomorphisms. 
Convergence degree, combinatorial and topological equivalence are introduced along with the 
ubiquitous Riemann surfaces known as jackets. Depending on our preference, we can utilize all 
or just a subset of the jackets to capture the degree, although in the most minimal case, one must 
sometimes introduce other 2-dimensional objects known as patches for completeness. Despite 
all this, we have just touched on the wealth of insights to be garnered from higher-dimensional 
cellular topology. 

The model discussed exclusively in the main text is the i.i.d. model for which the combinatorial 
and topological apparatus developed earlier was apt to characterize the amplitudes. Their 
scaling in the parameter N (the size of the tensor) is controlled by the degree of the underlying 
graph. Using the combinatorial and topological equivalence moves, we establish the concept 
of core graph at a given order. It becomes apparent that one can order two distinct 1/A^ 
expansions, combinatorial or topological, according to the core graphs of the appropriate type. 
The important result is the at leading order, only the spherical topology contributes. 

We made a brief excursion into the land of matrix models where we used the concept of 
factorization introduced earlier to partition the degrees of freedom of the tensor model into 
subsets attached to specific jackets (and possibly patches). These subsets take the form of matrix 
models embedded inside the tensor structure. This identification provides a staging post to 
subsequently apply matrix model methods directly to certain sub-sectors of a given tensor model. 
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Afterwards, we returned to analyze the leading order of our i.i.d. model. The leading order 
graphs all exhibited a characteristic 'melonic' structure and could be mapped to a class of 
trees. This re-indexing of the series allowed us to resum the dominant contributions. This 
was already enough to reveal the presence of a critical behavior leading to a continuum limit. 
With this result, who could not venture to speculate on its physical implications. Some insight 
comes from dynamical triangulations, which can give a more geometrical interpretation to graph 
configurations. The melonic sector displays similar behavior to the branched polymer phase. 

In addition, the leading order correlation functions also satisfy some Schwinger-Dyson equa- 
tions. These in turn form a Lie algebra another representation of which is carried by the space 
of rooted colored D-aiy trees. 

Before concluding, we developed a rather abstract ansatz to obtain classical solutions to the 
tensor models. It is based on the various possible factorizations of the {D + l)-colored graphs. 
Ultimately, the main aim is to reveal properties of the tensor models by expanding around 
non-trivial classical solutions. 

10.4 Towards future endeavours 

To bring us in a full circle, colored tensor models enumerate random walks where the walks 
are themselves simplicial pseudo-manifolds. The feature that instigated the breakthrough for 
tensor models was the coloring. It allows one to identify the topology of the graphs, enables 
al/N expansion ordered, at least roughly, by topological features. Moreover, it is the coloring 
once again that facilitates an in-depth analysis of the dominant behavior and a method to 
index the algebra associated to the quantum equations of motion. But where do we see the 
theory progressing from this point on. Clearly, developments will be divided into two broad 
categories: generic methods to analyze tensor theories and applications to pertinent physical 
and mathematical problems. 

For the first category, their are a large number avenues to pursue. Let us list a few. 

• We should refine our knowledge of the topological features encoded by the various graphs, 
e.g. via the definition of (colored) homotopy invariants. This has an number of side 
benefits. For the various extant tensor models, it would allow us to investigate the nature 
of their sub-dominant contributions. Moreover, it may allow us to develop models which 
weight graphs according to these topological properties. 

• At the moment, we possess a rather nice map between melonic graphs and a class of 
trees. Of course, sub-leading graphs do not possess this melonic property. However, 
we should examine whether they possess a defining characteristic that we might aid in 
our task of resumming these orders. In any case, upon resummation, we should like to 
check the position of the first critical point for these higher order contributions. A nice 
property of certain matrix models is that the sub-dominant contributions all diverge at 
the same critical point as the spherical sector [53]. This allows for a certain tuning of the 
parameters so that on passing to the continuum limit, one incorporates contributions from 
all 2d topologies. This is known as the double scaling limit. 

• The double scaling limit that we mentioned a moment ago is a remarkable achievement 
of the matrix model program. Since higher dimensional topology is so much richer in 
structure, generic tensor models are unlikely to possess just a single double scaling limit 
but rather double and multiple scaling limits. One possibility has already been investigated 
in [83]. 

• Another point of paramount importance is the analysis of the Schwinger-Dyson equations 
and Ward Identities. These not only encapsulate the quantum dynamics but succinctly de- 
scribe the quantum symmetries. At leading order in the 1/A^-expansion of the i.i.d. model. 
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these symmetry generators have been shown to form a Lie algebra. A pertinent question is 
whether this algebra forms a representation of the dif feomorphism algebra of the -D-sphere, 
or perhaps even a sub-algebra. 

• The wholescale incorporation of matrix model techniques into colored tensor models is 
certainly possible in the following sense. They may be utilized as-is for certain sub-classes 
of tensor models which weight only one particular Riemann surface embedded in each 
(D + l)-colored graph. One has enough control at the level of the action to ensure such 
a weighting. The task is then is to repeat the study done throughout this article for 
such models. They certainly have different contributions at leading order and therefore, 
different critical behaviors. While there are several technical issues pertaining to this 
analysis, one must also attempt to give these models some physical interpretation. 

On the other hand, various matrix model techniques are not obviously generalizable. For 
example, while there is certainly enough freedom in the choice of basis to diagonalize an 
Hermitian matrix, there is seemingly not enough freedom to diagonalize a tensor without 
first drastically reducing the number of degrees of freedom (essentially from A^-^ down 
to A^^). However, recent work in relating these tensors to D-ary algebras [122, 123, 125] 
might help to motivate the study of both issues. 

Of course, this is just a non-exhaustive list of active topics, but somehow the most immediate. 

The problems of the second category may be summarized in the manner of a recent (in) famous 
press statement: there are known knowns, known unknowns, unknown knowns and unknown 
unknowns. The later subcategories are truly the most populated. However, let us wax lyrical 
for a moment on some possibilities. 

• Given the obvious link between colored tensor models and simplicial (pseudo-)manifolds, 
a natural application centres on (simplicial) quantum gravity, or rather, simplicial quantum 
physics in general. 

With regards to quantum gravity, we described earlier recent endeavours to describe quan- 
tum space-time. The natural space-time interpretation of i.i.d. like models is via a map to 
dynamical triangulations, while the projector group field theory models are more directly 
interpreted via spin foams. The general thrust of colored tensor models differs slightly 
from that of group field theories although there is a high level of cross pollination. On 
the one hand, colored tensor models are well adapted to analysis via analogues of matrix 
model techniques, upon whose development this review has remained largely concentrated. 
On the other hand, group field theories, as models with genuinely infinitely many degrees 
of freedom from the outset, set their sights on incorporating techniques from standard 
quantum field theory, e.g. renormalization and so forth. 

Moreover, the incorporation of quantum matter is done slightly differently in colored 
tensor models and group field theories. For the former, the most readily feasible matter 
coupling is of Ising/Potts or hard dimer type [32, 34]. Once again, much inspiration 
can be taken from previous matrix models studies. For the latter, there is a tentative yet 
ambitious proposal that matter might be incorporated via topological singularities charged 
with the relevant quantum numbers of the matter in question [63, 105]. These topological 
singularities are naturally present in the Feynman diagrams of the theory - they put 
the pseudo into pseudo-manifold. This has been concretely studied in 3 dimensions with 
favorable results. We should investigate this idea in more physical regimes. 

• Matrix models have become a rather convenient tool with which to study certain problems 
in mathematical biology, in particular, the secondary and tertiary structures of ribonucleic 
acid (RNA) [49, 114]. RNA is a sequence of molecules, called nucleotides, which encodes 
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genetic information. The sequence itself is known as the primary structure. However, the 
molecules in the sequence are attracted to each other, which causes the sequence to fold. 
This folding is the secondary structure. Moreover, these folds arrange themselves not just 
on the plane but in 3-dimensional space, generating the tertiary structure. Matrix models 
take advantage of their built-in topological expansion to model the tertiary structure. 
Colored tensor models provide a natural 3-dimensional setting within in which to embed 
the various folded sequences and may give a more realistic picture for the structures in 
question. 

• As we motivated at the beginning, in our line of work it is often important to be able to 
count, and moreover, to count precisely. We saw in Section 7 that the leading order of 
a class of tensor models could be mapped to the space of colored rooted {D + l)-ary trees, 
which in turn aided the leading order resummation. No doubt as one investigates higher 
and higher orders, one shall encounter maps between different spaces of graphs. Thus, one 
can expect a flow of information among tensor models, graph theory and combinatorics. 

As one can see, the general formalism of colored tensor models is broad enough to be appli- 
cable in principle to a wide variety of challenges. The task remains, however, to pinpoint more 
exactly such situations. 
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